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ABSTRACT 

An analysis of the first set of low-redshift (z<0.08) Type la supernovae monitored by the 
Carnegie Supernova Project between 2004 and 2006 is presented. The data consist of 
well-sampled, high-precision optical (ugriBV) and near-infrared (NIR; Y JHK S ) light 
curves in a well-understood photometric system. Methods are described for deriving 
light-curve parameters, and for building template light curves which are used to fit 
Type la supernova data in the ugriBVY J H bands. The intrinsic colors at maximum 
light are calibrated using a subsample of supernovae assumed to have suffered little or 
no reddening, enabling color excesses to be estimated for the full sample. The optical- 
NIR color excesses allow the properties of the reddening law in the host galaxies to be 
studied. A low average value of the total-to-selective absorption coefficient, Ry ~ 1.7, 
is derived when using the entire sample of supernovae. However, when the two highly 
reddened supernovae (SN 2005A and SN 2006X) in the sample are excluded, a value 
Ry ~ 3.2 is obtained, similar to the standard value for the Galaxy. The red colors of 
these two events are well matched by a model where multiple scattering of photons by 
circumstellar dust steepens the effective extinction law. The absolute peak magnitudes 
of the supernovae are studied in all bands using a two-parameter linear fit to the decline 
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rates and the colors at maximum light, or alternatively, the color excesses. In both cases, 
similar results are obtained with dispersions in absolute magnitude of 0.12-0.16 mag, 
depending on the specific filter-color combination. In contrast to the results obtained 
from the comparison of the color excesses, these fits of absolute magnitude give Ry ~ 1— 
2 when the dispersion is minimized, even when the two highly reddened supernovae are 
excluded. This discrepancy suggests that, beyond the "normal" interstellar reddening 
produced in the host galaxies, there is an intrinsic dispersion in the colors of Type la 
supernovae which is correlated with luminosity but independent of the decline rate. 
Finally, a Hubble diagram for the best-observed subsample of supernovae is produced 
by combining the results of the fits of absolute magnitude versus decline rate and color 
excess for each filter. The resulting scatter of 0.12 mag appears to be limited by the 
peculiar velocities of the host galaxies as evidenced by the strong correlation between 
the distance-modulus residuals observed in the individual filters. The implication is 
that the actual precision of Type la supernovae distances is 3-4%. 

Subject headings: distance scale - dust, extinction - galaxies: distances and redshifts - 
galaxies: ISM - supernovae: general 



INTRODUCTION 



The discovery in the late-1990s from measurements of di stant Type la supernovae (SNe la) 



that the expansion of the Universe is currently accelerating (IRiess et al.1 Il998l : iPerlmutter et al 



19991 ) has given rise to an exciting new era in cosmology. The observations appear to require 
a previously unrecognized form of energy density ("dark energy") which today is the dominant 
constituent of the Universe. In rec ent years, observa t ions of SNe la to increasingly higher redshift s 
have confirmed this finding (e.g., IRiess et al.l 120071 : lAstier et al.l 120061 : IWood-Vasev et al.l 120071 ) . 
Measurements of the angular power spectrum of the cosmic microwave background radiation provide 
independent confirmation, favoring a flat geomet ry for the Universe with a dominant 73% of the 
energy density in the form of dark energy (e.g., see lSpergel et al.ll2007l ). Baryon acoustic oscillations 
(jEisenstein et al.ll2005l) and mea surements of the ratio of X-ray-emitting gas to total mass in galaxy 



clusters (Allen et al 



2004 . 120071 ) have also provided convincing confirmation of this picture. 



Effort is currently focused on determining the nature of dark energy by measuring an equation- 
of-state parameter of the form w = P/(pc 2 ), and its time derivative w. Several ongoing and 
future projects aim at filling the Hubble diagram with high-redshift (z > 0.1) SNe la in order to 
measure w and w with sufficient precision to distinguish among the various proposed models of 
dark energy. However, these measurements require an improved reference set of SN la observations 
at low redshift. Surprisingly, the currently available low-redshift datasets are neither sufficiently 
numerous nor homogeneous to adequately complement the high-redshift data. 



The Carnegie Supernova Project (CSP) is one of several ongoing efforts to improve the quality 
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of the low-redshift data0 Over the five-year period ending in May 2009, the CSP has obtained 
densely sampled optical and near-in frared (NIR) light c urves of ~100 SNe la in a well-understood, 
homogeneous photometric system (jHamuv et al.l 120061 . hereafter H06) . The NIR was included 
specifically to address the unknown intrinsic colors and interstellar reddening to SNe la. This work 
presents an analysis of the first release of SNe la by the CSP. 

In § [21 we present an analysis of the light curves, both optical and NIR, and introduce the 
methods used to derive parameters which serve to determine distances. In § [3j the non-trivial issue 
of measuring extinction is considered, and the nature of the reddening law in the host galaxies is 
studied. In § the properties of SNe la as standardizable candles is examined by using Hubble- 
flow distances to fit the relationship between absolute peak magnitudes, decline rates, and colors 
or reddenings for different bands. Finally, in § [U the implications of our findings are discussed. 

This paper is the second of three companion papers presenting first results from CSP obser- 
vations of low- and high-redshift SNe l a. The low-redshift d ata used in the present analysis are 
described in detail in the first paper bv lContreras et al.l (|2009i . hereafter C09). In the third paper, 
NIR light curves of 35 high-redshift SNe la observed with the Magellan Baade 6.5-m telescope 
are combined with the low-redshift data in order to construct the first rest-frame z-band Hubble 



diagram of SNe la (IFreedman et al.ll2009i ). 



2. PROPERTIES OF THE LIGHT CURVES 

The CSP carries out its follow-up observations with the Las Campanas Observatory (LCO) 
Swope 1-m telescope using optical ugriBV bands, and NIR YJH bands. Additionally, Y JHK S 
photometry is obtained with the Wide-field Infrared Camera (WIRC) at the du Pont 2.5 m telescope 
of LCO. As mentioned in H06, the NIR monitoring of the CSP during the first year of observations 
was generally not as thorough as that for the optical due to the fact that RetroCam, the NIR 
camera at the Swope telescope, was not available until September 2005. This situation was much 
improved during the final four campaigns of the project. 

The photometry of each SN is performed relativ e to nearby fie ld stars which are, i n turn , 
calibrated f rom observations of standard stars in the lLandoltl (|1992l ) (BV), ISmith et al.l ([2002 ) 
(ugri), and iPersson et al.l (|2004l ) (Y JHK S ) systems during several photometric nights. Point- 
spread function (PSF) photometry of the SN is measured after the underlying light from the host 
galaxy has been subtracted from every follow-up image. Details of the measurement procedures and 
photometric system are given by H06. Magnitudes are expressed in the natural photometric system 
of the Swope+CSP bands, as explained in C09. In accordance with that paper, in the remainder 
of this work we will refer to our optical natural photometry as ugriBV . The final photometry for 



9 See Contreras et al, ( 20091 ) for a summary of other groups that are producing large databases of low-redshift 
SN la light curves. 
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the SNe is published by C09. 

In the following, the methods employed for deriving light-curve parameters for the SN sample 
are presented, including the construction of template light curves for the ugriBVYJH bands. 



2.1. Direct Measurements 

The SNe la included in this analysis are listed in Table [TJ Twenty-six of these SNe have optical 
observations beginning before l?-band maximum light and with a cadence of a few days or less, which 
allows direct functional fits to their data. A subsample of nine SNe in this group also have enough 
YJH coverage to allow the same type of fits around maximum light. These fits provide direct 
measurements of light-curve parameters such as peak magnitude m max , time o f maximum t 



max 

and decline rate Amis(i?), the latter measured from the i?-band light curves |Philhpjll99d)E°l 
We therefore decided to directly measure the light-curve shape parameter, Ami^(B), whenever 
possible rather than derive it from template fits in order to tie our luminosity corrections to a 
purely empirical parameter. These resulting parameters are given in Table O labeled as "Spline" 
in column 7. The fits also serve to build a set of template ugriBVYJH light curves which may be 
used to fit more poorly observed SNe, as described in § 12.21 

Spline functions are preferred for these direct fits since they are more readily adaptable to 
the shapes occurring in S N la light cur ves than are other functions such as polynomials. We 



use the FITPACK library (IDierckxl Il993l ). which employs a variable number of knots depending 
on the sampling and precision of the data. These spline functions allow the whole time span of 
the observations to be fit without introducing undesirable oscillations where the coverage is less 
dense. The high quality of the photometry and the good time sampling of the observations ensures 
the accuracy of the fits. Errors in the parameters derived from the spline fits are modeled via a 
Monte Carlo calculation performed by randomly perturbing each data point according to its quoted 
photometric error and then measuring the scatter of the output measurements. 



As mentioned by lContreras et al.1 (J2009J), one can expect a certain degree of correlation among 
light-curve points introduced by the host-galaxy subtraction process because the same template 
image is used to subtract all of the follow-up images. Not taking the correlation into account 
may lead to underestimating the derived uncertainties of the light-curve parameters. The effect 
is expected to be more noticeable when the host-galaxy subtraction is more uncertain, and the 
photometric precision and coverage are larger. We have estimated this effect by computing the 
covariance matrix for the data points and including this in the Monte Carlo calculations. The 
covariance matrix is estimated by examining the photometry of artificial stars which are added to 
the subtracted images at fixed locations around the SN and which are set to have the same flux 



10 Throughout this paper we use the superscript or subscript "max" to denote the quantity measured at the time 
of maximum light in the corresponding band. 
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as the SN at each epoch. In general, these experiments have yielded results in agreement with the 
case in which the correlations are neglected. Only in a few cases have the uncertainties in the peak 
magnitudes increased by more than a few thousandths of a magnitude. Therefore, we consider this 
effect to be negligible for the current analysis. 

K-corrections are applied to the observed magnitudes to convert th ese to the re s t fram e. For 
this purpose , we utilize a libra ry of SN la template spectra provided by iHsiao et al.l (|2007l ) in the 
optical, and lHsiao et al.l (|2009l ) in the NIR. The K-corrections are computed simultaneously with 
the light-curve fits. The template spectra corresponding to the epochs of the observations are 
first warped by a smooth function of wavelength in order to make the synthetic colors match the 
observed ones for all available bands. Since the template spectra cover a limited range of epochs 
from —15 to +85 days relative to B-band maximum, observations outside these limits are ignored. 

The accuracy of the K-corrections can be tested by comparing them with values computed 
from the 194 follow-up spectra obtained by the CSP of the sample SNe. The spectra typically 
cover the optical griBV bands and a range of epochs from —12 to over +100 days with respect 
to .B-band maximum light. In most cases, the agreement with the K-corrections derived from the 
template spectra is better than ±0.02 mag. Unfortunately, the lack of spectral coverage in the u 
and NIR bands prevents this test from being performed in those bands. 

Figure [U shows the spline fits for all SNe with pre-maximum data in each band. It should be 
noted that since normal SNe la peak earlier in u and i than they peak in B (see § 12. 41 and Figure U]), 
some SNe were not observed before maximum light in these two bands, reducing the sample from 
26 in B, to 21 and 19 (respectively) in u and i. 

The value of Am\^{B) obtained for each SN is indicated next to the corresponding curves in 
Figure[H The range of decline rates covered by the sample is 0.7 < Ami^(B) < 1.9 for the BVgr 
bands, and 0.8 < Am\^{B) < 1.9 for the uiYJH bands, which corresponds well to the full range 
of observed values of Am\^{B) for SNe la. 

The optical data for these 26 SNe start between —12 and —1 days relative to B-band maximum 
light, and end between 17 and over 100 days past maximum, typically extending more than one 
month after maximum light. The YJH data for the nine SNe with pre-maximum coverage start 
between —10 and —1 days, and end between 5 and 90 days past maximum. 



Following the work of lHamuv et al.l (|1996bl ) and lElias- Rosa et al.l (|2008l ). we searched for alter- 
native characterizations of the SN la decline rates by comparing Ami^(B) with values of Arrit(X) 
for different bands X, and times t, since maximum. Figure [2] shows the cases for ugriBV which 
provide the strongest correlations with Ami5(-B); Figure [3] shows the same results for the NIR 
YJH bands. In general, at longer wavelengths, the best correlations with Ami^(B) are found at 
later epochs. Even in the case of the NIR iYJH bands, where the double-peaked shape of the light 
curves complicates matters, reasonable correlations with Am\^{B) are obtained at t ~ 35-40 days 
— i.e., during the radiative decline. 
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The behavior of the timing of maximum light for different bands with respect to B band was 
also studied. This is shown in Figure 01 For ugVr, the delay (or advance for it) is roughly constant 
for SNe of all decline rates. In the iYJH bands, there is a nearly constant advance of 2-4 days for 
SNe with Amis(5) < 1.5, whereas there is a delay of several days for SNe with /\m\^{B) rs 1.8. 
Note that the scatter in Figured] is generally greater than the error bars, implying that there is an 
intrinsic component. 



2.2. Template Light-Curve Fits 

The availability of such a complete set of ugriBVY J H light curves with pre-maximum cov- 
erage allows template light curves to be built which can be used to derive light-curve parameters 
for SNe with inadequate coverage around maximum brightn ess. For this pur pose, we used the 
SNOOPy (SuperNovae in Object-Oriented Python) package (IBurns et al.ll2009l), whi ch is derived 



from the the multi-filter light-curve fitting method developed by iPrieto et al.l (|2006) , with a few 
modificat ions. The techniqu e for generating templates is briefly summarized below, and in more 
detail bv burns et all (|2009l ). 



The SNe that were fit with spline functions (see § 12. ip form the training dataset. The time of 
£?-band maximum light t max (B), the decline rate measured as Ami^(B), and the peak magnitude 
m max ^ n eac h band X are used to place the photometric data points of the training set in a three- 
dimensional (3-D) space for each band defined by the following variables: 

1. The rest-frame epoch with respect to 5-band maximum, At = [t — i max (l?)]/(l + z), where t 
is the time of the observation and z is the heliocentric redshift of the SN; 

2. the decline-rate parameter, Am\^(B), considered as a single parameter to characterize all 
bands; and 



3. the K-corrected magnitude relative to maximum light, Am = mx(t) 
each observed magnitude mx{t)- 



in 



x 



K x (t), for 



SN 2006X was removed from the training dataset in ugriBV because of its peculiar behavior, 
especially in the B band, which showed an abnormally flat decline rate beginning one month after 
maximum light. This morphology has recently been expla ined as a light echo produced by dust 
surrounding this heavily reddened SN (jWang et al.1 l2008bl ). However, since there is no evidence 
that the NIR light curves of SN 2006X suffered such contamination, the data for this event were 
maintained in the training set for building YJH templates. 

The data points of the training set define a surface in the 3-D space (At, Am\^(B), Am) 
for each band. To produce a smooth interpolation to any point in the 3-D space, we use a 
2-D variation of the "gloess" algorithm, which is a Gaus s ian-w indowed and error-weighted ex- 
tension of data-smoothing methods outlined by IClevelandl ()1979l ) that was first implemented by 
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Persson et al.l (l2004r) fo r fittin g Cepheid light curves. This is th e main difference between the 



method of IPrieto et all (|2009 ) and that of the SNOOPy package (jBurns et al J 12003 ). While the 
former method uses the fit spline functions as templates and interpolates among them, SNOOPy 
performs an interpolation among the data points of the SNe. The gloess algorithm ensures a smooth 
interpolation even with heterogeneously sampled data. 

The SNOOPy package generates templates for a given value of Am\^{B) by making a slice 
through the (At, Arrii^(B), Am) surface, interpolating along a constant Am\§{B) line. The 
best-fit time of maximum light i max , decline rate Am\^{B), and peak magnitudes m ma,x are then 
determined by fitting these templates to the ugriBVYJH data of each SN via x 2 minimization, 
with K-corrections computed in an iterative fashion, as described in § 12.11 Uncertainties in the 
fit parameters are taken from the diagonals of the covar iance matrix outp ut from the Levenbert- 
Marquardt least-square routine used to fit the templates (jBurns et al.ll2009l ). Additionally, we have 
tested the possible effect of correlated error s in the photometry d ue to the use of a single host- 
galaxy image in the subtraction process (see IContreras et al.ll2009l ). Similarly to the results of the 
tests for spline fits that are given in § 12.11 in the case of template fits we find that the effect is very 
small in all but a few examples. We have thus decided to use the uncertainties that are derived 
from uncorrelated errors. 

This way, it is possible to derive the light-curve parameters for each SN, including the objects 
in our sample which were not observed until after maximum light. Table [2] lists the results of these 
template fits for the cases in which a spline function could not be fit; these are labeled as "Tempi." 
in column 7. The template fits themselves are shown by C09 together with the observed light curves 
for all SNe in the sample. 

In order to test the accuracy of the method, template fits were performed for the subset of well- 
observed SNe and the resulting parameters were compared with those obtained directly from spline 
fits (§ 12. ip . These comparisons are shown in Figure[5]for the peak magnitudes in ugriBVYJH, and 
for the decline-rate parameter Am\§{B). The average difference in peak magnitudes is ±0.05 mag 
or less for all of the optical filters except the i band, with no obvious systematic differences except, 
perhaps, in u. The dispersion in the iYJH bands is ~ ±0.1 mag, with some evidence for systematic 
differences in YJH. A similar situation holds for the differences in the time of maximum as 
measured with spline and template fits, where good consistency (±0.6 days or less) is found in 
ugrBV, but the iYJH bands yield worse agreement. The relatively poor precision of the YJH 
template fits is at least partly due to the small sample used to derive the templates. However, 
the large dispersion in i is the product of variations in the morphology of the secondary maximum 
(see § 12. 3h . Figure [5] shows that, except at the fastest decline rates, there is reasonable consistency 
between the decline rate as measured by spline and template fits. 

The uncertainties associated with extrapolating the templates to obtain peak magnitudes and 
decline rates for those SNe caught after maximum light are more difficult to quantify. Clearly, the 
later the observations begin, the larger are the possible extrapolation errors. One way to approach 
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this issue is to take the spline fits of the light curves of the SNe in the sample that were well- 
observed before and after maximum, resample these to simulate a SN whose observations did not 
begin until after maximum, and then fit these data using the ugriBVY J H templates to recover 
the time of B maximum, Amis(-B), and the maxim um-light magnitud es. Calculations of this type 
will be presented in more detail in a future paper ([Burns et al.l 120091 ) . but for the representative 
case of a SN for which coverage began ~1 week after maximum and continued until ~70 days after 
maximum, we find a typical random uncertainty of ~0.1 mag in the peak magnitude, averaged 
over all filters, and a much smaller systematic difference of ~0.03 mag. Experience shows that for 
SNe with photo metry beginning l ater t han this, template fits are not reliable. Fortunately, only 
one event in the IContreras et al.l (|2009l ) sample, SN 2004dt, falls into this category; hence, it is 
excluded from the analysis in the present paper. 

The lack of a sufficiently large set of SNe with good coverage around maximum light in the K s 
band prevents the production of template light curves as was done for al l other bands. Instead, we 
use the polynomial tem plate light curve introduced by iKrisciunas et al.1 (|2004al ). A stretch factor 
(jPerlmutter et al.lll997l ) based on the decline rate of each SN is applied to the time axis of the 
polynomial — along with the redshift-dependent time dilation — which is then fit to the data in 
order to obtain a peak magnitude. This kind of fit may be utilized for SNe with data in the range 
of validity of the template, which is between —12 and +10 days relative to t ma,x (B) for a SN of 
stretch factor S = 1. Th e specific formula employed to convert between Amis (5) and S is that 



given by I Jha et al.l (|2006l ): 



S 



3.06 - Am 15 (B) 
2M ' 



(1) 



2.3. The Secondary Maximum 



Hdflich et al.l (|1995l ) and lPinto &: Eastman! (|2000l ) were the first to attempt a t heoretical un 



derstanding of the NIR secondary maximum in SNe la. More recently, iKasenl (|2006l ) modeled the 
NIR light curves of SNe la, arguing that the timing and strength of the secondary maximum is a 
consequence of the ionization evolution of the iron-peak elements in the ejecta. A specific prediction 
is that more luminous SNe should have a later and more prominent secondary maximum. In order 
to test this, we examined the delay and relative strength of the secondary maxi mum in thej-band 
light curves of the CSP SNe with the best coverage. Using the nomenclature of IKasenl (120061 ). the 
following parameters were measured: 



1. The phase of the peak of the secondary maximum with respect to £>-band maximum, ti — 1&\ 

2. the difference in magnitudes between the secondary maximum and the local minimum between 
the primary and the secondary maxima, m,2 — mo; and 



3. the difference in magnitudes between the primary and secondary maxima, m\ — m%. 
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Figure [6] shows the results plotted as a function of Amis (5). If the decline-rate parameter 
Amifj(-B) (or peak luminosity) is associated with the amount of 56 Ni mass sy nthesized in the ejecta 



(|Arnettlll982l : lArnett et al.lll985l : iNugent et al.lll995l : IStritzinger et al.ll2006l ). then Figure © can be 
compared directly with Figure 11 of iKasenl ()2006l ). Indeed, the upper panel of Figure [6] shows 



a strong correlatio n between Am^(B) and t he timing of the secon dary maximum, confirming 
previous results by lHamuy et al.l ()1996bl ) and lElias-Rosa et al.1 (J2008J). However, the middle and 
lower panels of Figure [6] show that there is little evidence for a dependence between Amis (5) an d 
the strength of the secondary maximum, as mea sured w i th res pect to either the primary maximum 
or the local minimum between the two maxima. IKasenl (|2006l ) points out that both the timing and 
strength of the secondary maximum can also be affected by the outward mixing of 56 Ni in the ejecta, 
the amount of stable iron group elements produced in the explosion, the progenitor metallicity, and 
the abundance of calcium in the ejecta. Of these, mixing appears to have the greatest effect on 
the strength of the secondary maximum, and therefore may be responsible for the lack of obvious 
correlations in the middle and bottom panels of Figure [6l 



Krisciunas et al.l (j200ll ) suggested a different way of measuring the strength of the /-band 
second maximum which consists of converting the light-curve observations to fluxes normalized to 
the flux at the primary maximum, fitting a high-order polynomial to these, and then determining 
the mean flux (based on an integration of the polynomial) from 20 to 40 days after the time of B 
maximum. They found a good correlation between this parameter, (I)20-40j an d the decline rate. 
Similar measurements were carried out for the i-band light curves of the best-observed subsample 
of CSP SNe, and plotted vs. Amis (2?) i n Figure UJ Illustrated as a solid line are the same 
measurements for the family of SNOOPy templates. There is a clear trend in this diagram in the 
sense that luminous, slowly declining SNe generally display stronger secondary maxima. However, 
there is also a significant and real dispersion which is illustrated by the four labelled SNe — 2006D, 
2006bh, 2004ef, and 2004eo. The left half of Figure [8] shows that all four have very similar B and V 
light curves covering a narrow range of decline rates (Amis(-B) = 1.37-1.42). Nevertheless, there 
are significant differences in the strength and morphology of the i-band se condary max imum which, 
again, may reflect varying amounts of mixing of the 56 Ni into the ejecta ( Kasenl 12000 ) . 



These variations in the strength of the secondary maximum represent a significant problem for 
template fitting in the iYJHK bands, where the strength of the secondary maximum is assumed to 
be a smoothly varying function of Amis(-B). This is illustrated in the right half of Figure where 
the difference in the maximum-light magnitudes in BVi as derived from template and spline fits 
is plotted for the same four SNe (2006D, 2006bh, 2004ef, and 2004eo) with nearly identical decline 
rates. In the B band, the template fits give results that are in excellent agreement (±0.02 mag) 
with the direct measurements from the spline fits. The agreement in V is nearly as good, but 
in the % band, the template fits give errors as large as ~ ±0.15 mag. Figure [8] shows that these 
errors are a function of the strength of the secondary maximum: template fits to SNe with a weak 
secondary maximum will yield maximum-light magnitudes that are too faint, whereas fits to SNe 
with a strong secondary maximum will give maximum-light magnitudes that are too bright. This 
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effect most likely explains the relatively large dispersion in the difference in the peak i magnitudes 
measured from template and spline fits observed in Figure Until an observational parameter can 
be found which makes these differences predictable, they will remain a significant impediment to 
the fitting of template light curves in the iYJHK bands. Quantities such as (1)20-40 or ni2 — mo 
may prove to be such a parameter, but a larger sample of SNe is needed to test their utility. 



The results of the light-curve fits for the SNe in our sample are summarized in Table [2J where 
the times of maximum (i max ) and apparent peak magnitudes (m max ) are listed for all bands. Also 
given are decline-rate parameters (Am*) for each filter expressed as the amount in magnitudes that 
the light curve declines between maximum light and an epoch t days after maximum. Finally, the 
number of data points and range of epochs covered with respect to the time of maximum, as well 
as the fitting method (spline functions, § 12.11 or templates, § I2.2[) . are specified for each band. 

Based on the quality of the data and fits, we define a subsample of SNe with the best-measured 
light-curve parameters, mostly via direct fits but also through template fitting. Twenty-nine SNe 
belong to this group of "best-observed" events. These objects are identified in Table [TJ in Table [21 
their names are set in bold face. Special use will be made of this subsample when studying the 
relation between peak luminosity and decline rates in § |U 



from dust maps of the Galaxy may be used to correct for the interstellar reddening produced in our 
Galaxy. Beyond the Milky Way, the situation is considerably more complicated as there are at least 
three separate potential sources of reddening due to (a) dust in the intergalactic medium, (b) dust in 
the interstellar medium of the host galaxy, and (c) dust associated with the circumstellar material 
of the SN progenitor. We have no a priori knowledge of the dust properties in these different 
environments and the reddening laws could, in principle, be different for all three. Nevertheless, by 
comparing the colors of SNe with similar decline rates and taking advantage of the wide wavelength 
coverage of our observations, one can hope to discern, in general terms, the combined effect of these 
different sources of reddening, which we lump together and refer to as the "host-galaxy reddening." 



2.4. Light-Curve Parameters 



3. 



HOST-GALAXY REDDENING 



The E(B — V)cai values listed in Table CD derived 




3.1. Low-Reddening Sample 



In order to use color informatio n to determine redde ning of the SN light produced in the host 
galaxy, we follow the procedure of Phillips et al. ( 19991 ) of selecting a subsample of SNe which 
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are suspected to have suffered little or no dust extinction. This subsample will be used to derive 
intrinsic colors for SNe la as a function of decline rate which, in turn, will serve as a reference for 
computing color excesses for the whole SN sample. 

The criteria employed to select the low-reddening subsample of SNe are (a) SNe which occurred 
in E/SO-type galaxies, or SNe located away from the arms or nuclei of spiral galaxies, and (b) 
absence of detectable interstellar Na I D lines in early-time spectra. 

Ten SNe from the sample fulfill both criteria: SNe 2004eo, 2005M, 2005al, 2005am, 2005el, 
2005hc, 2005iq, 2005ke, 2005ki, and 2006bh. In what follows, these objects are used to derive in- 
trinsic color relations as a function of the decline-rate parameter Ami^(B). We use colors corrected 
for Galactic reddening since such correction is, in most cases, small (see Table [1]) and known to 
higher accuracy than the host-galaxy reddening. We thus do not perform any cut based on Galactic 
reddening. We assume a value of Ry = 3.1, which may not be exactly correct for each SN. How- 
ever, since the maximum Galactic reddening for this sample of ten SNe is E(B — V) = 0.14 mag, 
uncertainties in this correction will be small. 



3.2. The B - V Tail 



The subsample of low-reddening SNe can be used to study the behavior o f the B — V color 
approximately one month after maximum light. It was first found by iLiral ()1995l ) that unreddened 
SNe la follow a very similar linear regime in B — V between one and three months after maximum. 
As shown in Figure O we confirm this result and derive the following intrinsic color law for the 
"tail" of the B — V color evolution: 



(B - V) = 0.732(0.006) - 0.0095(0.0005) (ty - 55). 



(2) 



For consistency with Phillips et al.l dl999j), rest-frame days since maximum V light, ty, is used for 
the time axis. Equation (|2|) is valid over the range 30 < ty < 80 days and yields a dispersion of 
0.077 mag. The SNe employed in the fit cove r a broad rang e of de cline rates, 0.85 < Amis (5) < 
1.76 mag. This relation is similar to the one of lPhillips et al.l (|1999l ). who found, for a sample of six 
SNe and a different photometric system, (B — V)q = 0.784—0.0118 (ty — 55) with a dispersion of 0.06 
mag. Figure [9] shows that the residuals for individual SNe are strongly correlated, which implies 
that some of the 10 SNe in this subsample actually do suffer measurable host-galaxy reddening, or 
that there is an intrinsic dispersion in the Lira relationship, or both. We also note that the three 
SNe in our sample which occurred in early- type galaxies (E/S0), namely SNe 2005al, 2005el, and 
2005ki, show remarkably similar B — V color evolution — and also in other optical colors — during 
the tail, with an average deviation of —0.05 ± 0.003 mag from equation ([2]), and a dispersion of 
~ 0.035 mag. 



Equation ([2]) can be used to estimate color excesses, E(B— V)taii> for the SNe in our sample with 



-12- 



observations in the tail0 The results are listed in column (2) of Tabled! The uncertainties in these 
color excesses reflect only the photometric errors in the data points. In considering the color excess 
for any particular SN, the observed dispersion of 0.077 mag in the Lira law should also be added 
in quadrature to these errors. However, in the analysis presented in § 14.21 of the precision to which 
SNe la may be used as standardizable candles, this dispersion is included in the cjsn values that are 
yielded by the fits, and the errors used for the color excesses are those given in Table HI Note that 
since equation ([2]) was derived by fitting to the average color evolu tion of the low- r edden ing sample, 
it is guaranteed to generate some small negative color excesses. iPhillips et al.l (|1999l ) dealt with 
this by applying a Bayesian prior consisting of a one-si ded Gaussian distribution of A r values with 
a maximum at zero and a = 0.3 mag. More recently, iJha. Riess. &: Kirshnerl (|2007l ) used a prior 
consisting of a Gaussian with cjib-v) = 0.068 mag, representing the intrinsic dispersion in color, 
convolved with an exponential reddening distribution with a scale length t e ^ b _ v ^ = 0.138 mag 
which they derived from fits to nearby SN la light curves by assuming a particular intrinsic color 
distribution. However, until it is possible to confidently separate dust reddening from intrinsic color 
variation, we choose not to apply a prior to the color-excess measurements. 

For each SN, we tested the agreement of the B — V tail slope with that of equation ([2]). Good 
consistency was generally found, with the notable exception of SN 2006X, the most highly reddened 
SN in the sample. As mentioned in § 12.11 this SN showed a peculiar beh avior in the B band at the 
epochs considered here, probably due to the appearance of a light echo (jWang et al.ll2008bl ). 



Figure [TO] displays a plot of E(B — V)ta,n versus the equivalent width of the absorption in 
the Na I D lines produced in the host galaxies as measured from spectra we obtained of the SNe. 
Generally speaking, these spectra are of low wavelength resolution, and the Na I D lines are only 
visible when the equivalent width reaches ~1 A or more. The red symbols with dashed error bars in 
Figure [T0l indicate the low-reddening sample. Four of the eight SNe with the largest color excesses 
show detectable Na I D absorption, and the two with the largest color excesses, 2004gu and 20051u, 
correspond to the two with the highest Na I D equivalent widths. However, two SNe with slightly 
negative values of E{B — V)x,ax\, SNe 2005bg and 2006mr, also show strong Na I D absorption. In 
the case of SN 2005bg, the estimate of E{B — V)t a ii is based on a single B — V color measurement 
in the tail, and therefore may be unreliable. SN 2006mr was a very rapidly decli n ing e vent which 
occurred in the dusty central regions of Fornax A (NGC 1316). iGarnavich et al.1 (120041 ) studied a 
small sample of fast-declining SNe la and found that they followed the Lira relationship fairly well, 
but SN 2006mr suggests that there may be exceptions. In an y case, the correlati on between color 
excess and the Na I D equivalent width is known to be poor (jBlondin et al.l l2009). 



Figure OH shows that five additional SNe (2004ey, 2005ag, 2006D, 2006ax, and 2006gt) have 
color excesses E{B — V\ a n < 0.08 mag and undetectable Na I D absorption. In the following 
analysis, these five SNe are added to the original subsample of ten assumed to have suffered little 



11 We remind the reader that equation ([2]) and all other relationships given in this paper which involve the CSP 
observed magnitudes are valid only for photometry in the natural system of the Swope+CSP bands. 
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or no host-galaxy reddening. 



3.3. Pseudocolors at Maximum Light 

Since the CSP data — and generally most SN la observations — have better coverage around 
maximum light than during the declining tail, it is important to develop methods to determine color 
excesses from observations around the time of maximum. Using the peak magnitudes (corrected 
for Galactic reddening and the K-correction) derived in § pseudocolor^! can be computed as 
the differences between the maximum- light magnitudes in two given bands. The dependence of 
these pseudocolors on the decline rate, parameterized by Amis (B), can then be investigated for a 
subsample of SNe with low reddening. 

Figure [TT] shows the pseudocolors as a function of Amis (2?) for several band combinations 
which will be used in the analysis of § [4] The solid points in the figure correspond to the 15 SNe 
with low reddening. As expected, in all colors these SNe lie at the lowest values, following a roughly 
linear trend with decline rate over the range 0.8 < Amis (2?) < 1.7 mag. At the fastest decline 
rates (Amis (2?) > 1.7 mag), however, the pseudocolors at maximum light are significantly redder 
than these linear trends would predict. Table [3] summarizes straightdine fits to the colors of the 
low-reddening SNe with decline rates 0.8 < Amis (23) < 1.7 mag; these fits are plotted as dashed 
lines in Figure [TT1 Note that the slopes of the relations are largest when the uYJH bands are 
involved. In general, there is a correlation between the dispersion in color and the difference in 
the effective wavelengths, AA, of the two filters that form the color. Within this trend, the colors 
involving the u filter have a relatively higher dispersion for their AA values, whereas the colors 
involving the J filter have a lower dispersion. 

In the case of the (2? max — y max ) pseudocolor, the following linear law for the 14 SNe with low 
reddening and 0.8 < Amis (2?) < 1-7 mag is found: 



(B 



max t /max 



V 1 



-0.016(0.014) + 0.12(0.05) [Am 15 (B) - 1.1]. 



(3) 



Th e root-mean square (rms) of the fit is 0-0 6 mag . This relation is slightly redder than that given 
bv lPhillips et al.1 (|1999l ) and lAltavilla et all (J2004J). We use this intrinsic-color law to derive color 
excesses E(B — V) max for all SNe in the range 0.7 < Amis (2?) < 1-7 mag, computed as (2? max — 
ym&x^ _ ^max _ ymaxj^ rj,^ resu ^j n g values are listed in column 3 of Table [4] The uncertainties 
in these color excesses are calculated from the photometric errors in the peak magnitudes. 



The top-left panel of Figure [12] shows a comparison of the E(B — V) max measurements with 



The term pseudocolor is employed to stress the fact that these quantities do not represent the actual color of the 
SNe at any time but are just the difference between the magnitudes at maximum light of two bands, which occur at 
different times. 
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the E(B — V)taii values from § 13.21 The sample of 21 objects consists of the best-observed SNe 
for which both color excesses could be measured. In addition, the highly reddened SN 2006X has 
been excluded since it clearly deviates from the behavior of the rest of the SNe in the B — V tail, 
as mentioned in § 13.21 One can see that there is a small but systematic difference between these 
two color-excess measurements which, unexpectedly, appears to be correlated with the magnitude 
of the color excess. A formal fit to the data (plotted as a solid line in Figure fl~2l) gives E(B — 
V^taii - E ( B - ^)max = 0.032 - 0.496 - y) max , with a dispersion of 0.06 mag. The top-right 
panel of Figure [T2l shows that this difference also seems to correlate with the color of the SN. Since 
E(B — V)ta,ii and E(B — V) max are "observed" color excesses, as opposed to "true" color excesses 



(see iPhillips et al.lll999l ). and are measured at different epochs in the color evolution of the SN, 
one would expect there to be a small systematic difference between the two that would correlate 
with the SN color. This effect is produced primarily because the effective wavelengths of the filters 
evolve due to the significant color evolution of the SN spectral energy d istribution, but the re is also 



a smaller dependence on the total amount of dust reddening itself. The lHsiao et al.l (|2007l ) spectral 
template may be used to calculate the expected magnitude of these combined effects, which is 
plotted as a dotted line in both of the top panels of Figure [12] The observed systematic differences 
are greater than can be explained in this way, although the significance is only at the ~ 2a level. 

Since the method of estima ting the intrinsic co lors at maximum by bootstrapping from the 



Lira relation was first utilized bv lPhillips et al.l (|l999i ), we have examined the color excesses derived 



by these authors. Interestingly, the same trends observed in the top two panels of Figure [12] are 
clearly present in those measurements. Not only does this appear to confirm the reality of the 
effect, but it also argues that it is not an artifact of the CSP dataset. 

The lower-left panel of Figure [12] shows that E(B — V) ta ii — E{B — V) max does not depend 
significantly on the decline-rate parameter Ami^(B). The measurements suggest that there may 
be a weak correlation with the absolute B magnitude (corrected for decline rate), but this requires 
confirmation by more data. We do not currently understand the underlying causes of these corre- 
lations, but the fact that they exist implies that there is a component of one or the other of these 
color-excess measurements that is not due to interstellar dust. 

Keeping this caveat in mind, the fits corresponding to (y max — X™ ax ) for X\ = ugriBY JHK S 
can be used to derive color excesses E(V — AT^) max for all SNe with decline rates in the range 
0.8 < Ami5(.B) < 1.7 mag. Columns 4 through 11 of Tabled] list these color excesses. 



3.4. The Reddening Law 



The availability of photometric data spanning the u through K s bands enables the properties of 
the host-galaxy reddening t o be examined. Assuming that the dust extinction in the host galaxies 
obeys the law introduced by lCardelli. Clayton. &: Mathisl (|1989i ) and modified by lO'Donnelll (11994] ) 
(hereafter referred to as the "CCM+O law" ), we can investigate which value of the total-to-selective- 



-15- 



absorption coefficient Ry is favored by SN la data. Extinct ion in the Galaxy is well a pproximated 
by the CCM+O law with an average value of Ry ~ 3 (e.g., iFitzpatrick fe Massal l2007). 



3.4-1- Optical-NIR Color Excesses 

An average value of Ry for the subsample of best-observed SNe can be derived by comparing the 
color excesses between the optical and NIR bands presented in § 13.31 with the values of E(B — V) max . 
The use of optical-NIR color excesses ensures the necessary leverage in wavelength to provide a 
more precise measurement of the reddening law, even for SNe with low reddening. In Figure [T3j 
comparisons of E(V — A) max with X = iYJH are shown. 

Fits to the slopes E(V -X) max / 'E(B-V) 

max were carried out for all four band s considered here. 



The s lope were then converted to values of Ry, using equations (A5) and (A6) in iKrisciunas et al 



(|2006l ) and the CCM+O law coefficients ax and bx derived as explained in Appendix [Bj Table [5] 
summarizes the results of these fits. The data favor a low value of Ry, with an average of Ry = 
1.69 ±0.05 from the four fits. The corresponding slope for this average is indicated by the solid red 
lines in Figure [T3l 

If, however, the two highly reddened SNe, SN 2005A and SN 2006X, are excluded from these 
fits, the uncertainties increase but the favored values of Ry become larger and compatible with the 
standard Ry sa 3. An average of Ry = 3.2 ± 0.4 is derived from the four fits, and is indicated by 
the dashed blue lines in Figure [T3l 

This is a very interesting result which may indicate that SNe with low or moderate reddening 
are affected by standard dust extinction, while very highly reddened SNe suffer relatively more 
reddening than extinction, as inferred from a low value of Ry. We will return to this finding after 
analyzing the calibration of SNe la as standardizable candles in § [H 



3.4-2. Highly Reddened Supernovae 

Highly reddened objects provide the most precise information about the behavior of extinction 
as a function of wavelength. Our sample includes two objects, SNe 2005A and 2006X, for which 
E(B — V) > 1.0 mag. Although the optical light curves of SN 2 006X showed eviden ce of a light 
echo, this did not appear until ~1 month after maximum light (jWang et al.l l2008bl ). and so the 
peak magnitudes should be uncontaminated. Figure HH shows the values of E(V — X\) derived for 
both SNe (see § I3.3[) plotted as a function of the effective wavelength of each filter X. Note that at 
these extreme values of the color excesses, it is necessary to conv ert the observed valu es to true color 
excesses, and to also correct Ami^(B) for reddening effects (see iPhillips et al.lll999l ). Table [6] gives 
the linear relation s for converting from observed to true color excesses that are derived using the 
Hsiao et al.1 (|2007l ) template spectra. Using the same spectra, the following relation for correcting 
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the observed decline rate for reddening is found: 



Ami 5 (B) t 



Ami 5 (5) observed + 0.065E(B - V) observed- 



(4) 



Assuming that the reddening of SNe 2005A and 2006X is produced by dust that is described by 
the CCM+O reddening law, fits to the color excesses in Figure [T4"l give values of Ry of 1.68+0.10 and 
1.55 ± 0.07, respectively. These fits are plotted as solid lines in Figure and are to be compared 
with the dotted lines which show the best-fitting CCM+O reddening law where Ry is fixed at its 
canonical Galactic value of 3.1. Although the CCM+O law does a good job in th e optical and NIR, 
it fails miserably in fitting the (V — u) color excesses. Recently, iGoobarl (|2008l ) has considered a 
model where multiple scattering of photons by circumstellar dust steepens the effective extinction 
law, converting it to a power law. When this model is fit to the color excesses of SNe 2005A and 
2006X, excellent agreement at all wavelengths is obtained with Large Magellanic Cloud (LMC) 
dust and a power-law index p « —2.4. These fits are plotted as dashed lines in Figure [21 Thus, 
the data would appear to favor a model where the reddening of these two objects is produced by 
dust that is local to the SN. In this context, it is interesting to note that short-term variations in 



the in terstellar Na I D absorption lines in the spectrum of SN 2006X were observed by iPatat et al 



(j2007l ) and interpreted as evidence for the existence of significant circumstellar material. 



SN 2006X has also been studied in detail by IWang et al.1 (|2008al ). who found a best fit to 
the data with Ry = 1.48 ± 0.06, in excellent agreement with our results. However, these authors 
concluded that the d ata, includin g their [ /-band photometry, were matched quite well by the CCM 
reddening law. The IWang et al.l (|2008al ) color-excess measurements for SN 2006X were derived 
through comparison with a few individual SNe as well as via the method we employ here of looking 
at the maximum-light colors as a function of the decline rate. In general, the re is good ag r eement 
between their measurements and our own except for the critical U band, where IWang et al.l (|2008al ) 
find E(U — V)true = 2.65 ± 0.14 but we derive E(u — V)t rU e = 3.11 ± 0.14. Some of this discrepancy 
may be due to the fact that the sample of SNe la with well-observed light curves in the U or u 
bands is still fairly small, although in looking at Figure [TTJ it is difficult to see how the value of 
E(u — V) could be overestimated by 0.4 mag. Although the u filter response function is the least 
well-determined of the CSP bands, the color excesses for SNe 2005A and 2006X were derived with 
respect to other SNe la observed with exactly the same filter, and in the same photometric system. 
Thus, the color excess is well determined. The largest uncertainty is the effective wavelength that 
is associated with the CSP u band, which we estimate to be < 200 A. This maximum possible error 
is indicated by the horizontal error bars on the u points in Figure [T4l The conclusion is that the 
uncertainty in the u-band response function cannot explain the discrepancy. 

Only two other objects in our sample, SNe 2005kc and 2006eq, have significant host-galaxy 
reddenings and observ ations in the Y JH bands. Fits to the color excesses of these two SNe with 
both the CCM+O and lGoobarl (|2008l ) reddening laws are displayed in Figure [T5l The observations 
of SN 2005kc are consistent with either the CCM+O law with Ry = 4.4 ± 0.6, or a Goobar power 
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law with p = —0.7 ± 0.2. In the case of SN 2006eq, the measurements may be fit with either a 
CCM+0 law with Ry = 1.1 ± 0.8, or a Goobar power law with p = —3.5 ± 1.1. Thus, the results 
for SN 2006eq are similar to those obtained for SNe 2005A and 2006X, but the reddening law for 
SN 2005kc appears to be quite different, perhaps even consistent with normal Galactic reddening. 
These results confirm those of Figure [T3"l where SNe 2005kc and 2006eq are the two points with 
highest E(B — V) values among the moderately reddened SNe. These two SNe illustrate the power 
of optical-NIR color excesses in differentiating the reddening law. More well-observed, moderately 
reddened events are clearly needed to advance further in the study of the reddening laws of SNe la, 
and to look for possible correlations with other properties of the SN or its environment. 



TYPE IA SUPERNOVAE AS STAND ARDIZ ABLE CANDLES 



In this section, the homogeneous sample of CSP light curves is used to re-evaluate the precision 
to which SNe la may be used as standardizable candles. Two different approaches for calibrating the 
ab solute magni tudes of SNe la are considered. We start with the two-parameter method introduced 
by iTrippI (119981 ) in which the correlation between absolute peak magnitude, decline rate, and color 
is modeled for different choices of magnitudes and colors. We then move on to use the color excesses 
determined in § 13.4. II to examine the correlati on of reddening-corre cted, absolute peak magnitudes 
versus decline rate, following the approach of lPhillips et al.1 (|1999l ). The former method makes no 
assumptions about the intrinsic colors of SNe la; it simply considers decline rates and colors as two 
parameters that serve t o improve the pr e cision of SNe la as distance indicators. This formalism was 
originally employed bv lTripp &: Branch! (119991 ) to derive cos mological pa r amete rs from SN la data, 
and is part of the SALT (|Guv et al.l l2005h method used by lAstier et al.l (|2006r ) . The color term in 
this type of analysis makes no distinction between differences in intrinsic colors or in the amount 
of reddening suffered by the SNe. On the other hand, the second approach explicitly assumes that 
one can derive color excesses by determining the intrinsic colors of SNe la, and apply these as 
extinction corrections to the observed peak magnitudes. 

To calculate luminosity distances for the SNe, the following approximation is employed which 
is valid to a precision of better than 0.25% at redshifts z < 0.1 in a two-component (matter + dark 
energy) model of the Universe: 



4(zcmb; H , Om, ^a) 



(1 + ^helio) C 

(1 + ^cmb) Hq 



Z CMB + 7^ (^A Y + 1 ) ZCMLS 



(5) 



where Zhelio is the heliocentric redshift of the host galaxy, zqmb is the redshift referred to the 
CMB rest frame, and the cosmological p arameters were set to the standard values Hq = 72 km 
s- 1 Mpc" 1 , tt M = 0.28, and ft A = 0.72 (|Spergel et alJl2007l ). The factor (1 + z hc i io )/(l + z C mb) 
accounts for the fact that the photon redshift is observed with respect to the heliocentric reference 
system. The heliocentric redshifts for the sample SNe are given in Tabled! Also given in this table 
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are the redshifts in t he CMB frame calc ulated from the heliocentric redshifts using the velocity 



vector determined by iFixsen et al.l (|1996l ) from the COBE data (and which are available via the 
NASA/IPAC Extragalactic Database, NED). The uncertainty in the redshifts due to peculiar ve- 
locities is assumed to be a z = 0.001 (300 km s _1 in velocity). With the exception of the three SNe 
listed in Table [7] with z < 0.01 that were hosted by galaxies with direct distance measurements 
based either on Cepheids or the surface brightness fluctuation (SBF) method, we restrict the sample 
to those SNe with redshifts z > 0.01, where peculiar recession velocities become at least one order 
of magnitude smaller than the cosmological velocities. 

In the analyses presented here, the observed peak magnitudes of Table [2] are corrected for 
Galactic reddening using the values of E(B — T^Gai given in Table [TJ and adopting the CCM+O 
extinction law with Ry al = 3.1. 



4.1. Calibration Using Decline Rates and Colors 



The two-parameter model of iTrippI (jl998l ) assumes that the distance modulus of a SN has the 



following dependence on decline rate and color: 



fix = m x - M x {0) - b x [Am 15 (B) - 1.1] - p r x (my - m z ), (6) 

where mx is the observed peak magnitude in a given band X, and (my — mz) is a pseudocolor at 
maximum brightness from any choice of bands Y and Z. The fit parameters are the peak absolute 
magnitude of SNe la with Amis(S) = 1.1 and zero color, Mx(0); the slope of the luminosity vs. 
decline-rate relation, bx] and the slope of the luminosity-color relation 0x Z . Using measurements 
of peak magnitudes, pseudocolors, and decline rates {mxf, Ami5(B)f, (my — mz)i} (i = 1, • • • , N) 
for a sample of N SNe, the best-fit parameters are solved for via \ 2 minimization, as explained in 
Appendix |Al This is equivalent to minimizing the dispersion in the Hubble diagram. Note that a 
constant term has been added to the measurement uncertainties that appear in the denominator 
of the x 2 m equation (IA4h to account for possible intrinsic dispersion of the SN la data about the 
model given by equation Q. 



The lTrippI (|1998l ) model has traditionally been applied to fits of B versus Ami^(B) and B — V. 



Such an analysis is presented here, along with a fit of J versus Amis (5) an d V — J. 

Table [8] summarizes the results of the fits. Fit 1 was done using the entire sample of 32 SNe 
with reliable distances and the combination of B versus Ami^(B) and B — V. This fit is shown 
in panel (a) of Figure [TBI and yields a total dispersion of 0.17 mag and an intrinsic dispersion 
csn = 0.12 mag. Performing the fit with the subsample of 26 best-observed SNe (Fit 2) reduces 
the scatter to 0.15 mag and the derived intrinsic dispersion to just 0.09 mag without modifying 
substantially the fit parameters. This is shown in panel (b) of Figure [16j 
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Note that the fits described above include SNe with all the available decline rates; specifically, 
they include the fast-declining SNe with Am l5 (B) > 1.7 (SN 2005M, SN 2005ke, and SN 2006mr) 
which are marked with red circles in Figure [161 These SNe follow the behavior of the rest of the 
sample. This is substantiated by the results of Fits 3 and 4, done excluding these three SNe; 
the resulting fit parameters are in agreement within the uncertainties with those of Fits 1 and 2, 
and the dispersions found are nearly the same. This is a remarkable result since the red colors of 
the fast-declining SNe la — which occur preferentially in E and SO galaxies — are clearly largely 
intrinsic in origin, whereas the reddening of slower-declining SNe — which are most common in 
spiral galaxies — is presumably due at least in part to dust. 

If it is assumed for the moment that the pseudocolors (after correction for Galactic reddening) 
vary due to reddening by dust in the host galaxy, the color term, /?, can be converted to a total- 
to-selective absorption coefficient, Ry, adopting the CCM+O reddening law and following the 
prescriptions of Appendix [Bl The value of Ry obtained from each fit is given in column 7 of 
Table EJ As is seen, all the fits thus far considered yield values of Ry ~ 1.5, which is significantly 
different from the standard Galactic value of 3.1. In the right-hand panels of Figure [TBI the fits 
are labeled with the corresponding Ry values. For comparison, the reddening vector predicted by 
Ry = 3.1 is also plotted dotted line. 

The Tripp (1998) model was also fit to the subset of best-observed SNe that meet the condition 
that (S max — y max ) < 0,4 ma g. This excludes both the fast-declining objects, which are intrinsically 
red, and SNe 2005A and 2006X, which suffered heavy reddening due to dust. Fit 5 in Table [8] gives 
the results of such a fit. Once again, the zero point and slope of the luminosity vs. decline-decline 
rate relation (Mg(0) and &b) do not change significantly in comparison with the previous fits. The 
luminosity-color slope, /3g V , does change slightly (by ~ la), resulting in an even lower value of Ry. 
These results are relevant to the study of high-redshift SNe la for which red (and ther efore faint) 



objects [with (£? max — l/ mffit ) > 0.4 mag] are rarely observed (e.g.. iKowalski et al.ll2008l ). 



The results of applying the Tripp (1998) model to J versus Am\^{B) and V — J for the sample 
of best-observed SNe with J-band coverage (21 objects) are given in Fit 6 of Table EJ and are 
plotted in the lower part of Figure HH1 Two additional fits are presented in Table [8] for the J versus 
Amis (2?) and V — J combination. In Fit 7, the two fast-declining SNe with coverage in J are 
excluded. A slightly lower value of the slope py J is obtained, with an almost negligible intrinsic 
dispersion of ctsn = 0.02 mag. Fit 8 employs a color cut of (y max — j max ) < q q mag, which is 
equivalent to limiting the sample to (U max — y max ) < 0.4 ma g. Again, the best-fit parameters are 
in good agreement with those of Fit 7, and the dispersions are likewise very similar. 

Note that the interpretation of ctsn as the intrinsic dispersion of the SN data about the model 
relies on the correct estimation of the measurement uncertainties and their covariances. An error 
in the treatment of the uncertainties would lead to different usn- Consequently, the fit parameters 
themselves could potentially change. In order to test this, the dependence of the fit parameters on 
the adopted uncertainties in the peak magnitudes, colors, and decline rates was examined. This was 
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done by alternatively adding an arbitrary constant uncertainty on each quantity, up to 0.06 mag. 
While the values of <tsn obtained were found to decrease with increasing uncertainties as expected, 
the best-fit values of the parameters remained well within the uncertainties. 

In order to test the dependence of our results on the method used to derive light-curve parame- 
ters, we repeated the fits using peak magnitudes, decline rates and colors from template light-curve 
fits (see § I2.2p . We used the sample of best-observed SNe for which we performed both spline and 
template fits. The results are equivalent within the quoted errors to the ones presented above. 

4.2. Calibration Using Decline Rates and Extinction 

In this section, the color excesses derived in § [3] are used to produce fits of absolute peak 
magnitude versus decline rate and reddening. It is explicitly assumed that the color excesses, 
E(Y — Z), are due to dust and can be converted into an absorption in the X band via Ax = 
R x z E(Y — Z). This absorption is used to correct the observed peak magnitudes according to the 
model 

fix = m x - M x (0) - b x [Amis (B) - 1.1] - R Y X Z E{Y - Z). (7) 

Note the similarity of this model to that of equation ([6]) of § 14. 1\ especially if one considers that the 
color excesses derived in § 13.31 are based on linear relationships between the observed pseudocolors 
and Ami5(5)J^§ Since the color excesses were obtained for SNe in the range 0.7 < Amis(-B) < 1.7 
mag, the fast-declining SNe are necessarily excluded from the analysis of this section. 

We first focus on fits done with B and E{B — V) data. To provide a reference, the sample was 
initially limited to the twelve SNe which were previously identified as having suffered little or no 
dust extinction (see § I3.2p . These low-reddening objects were used to fit the model of equation ([7]) 
with the last term on the right-hand side set to zero — that is, without applying any reddening 
correction. The resulting parameters are given in Fit 1 of Table [H These twelve SNe with low 
reddening yield a total dispersion of 0.19 mag, with an intrinsic dispersion ctsn of the same amount. 

Next, the entire SN sample was fit by first fixing Ry to the standard Galactic value of 3.1. 
The results are given in Fit 2 of Table where we solve only for M#(0) and bs- This fit implies 
substantially higher peak luminosities and, most notably, a very large dispersion of ~ 0.5 mag. 
Thus, correcting for the observed reddening as if it were produced by typical Galactic interstellar 
dust yields a dispersion which is a factor of ~2.5 greater than that obtained for the subsample of 
low-reddening SNe. 

Beginning with Fit 3 of Table El R% v was allowed to be a free parameter. Here the full 



Indeed, as discussed in Appendix^ the two treatments are mathematically equivalent. 
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sample of SNe with 0.7 < Amis(5) < 1.7 mag was employed, while in Fit 4 only the best-observed 
subsample was considered. The resulting parameters for these two fits are indistinguishable within 
the errors. The best-fit value of Rg V in both cases is ~2.8, corresponding to Ry ~ 1-5. These two 
fits are shown in panels (a) and (b) of Figure [T71 

The values of Rq derived in Fits 3 and 4 are strongly influenced by the two highly reddened 
SNe 2005A and 2006X. Fits 5 and 6 in Table [9] give solutions where these two SNe are excluded 
from the sample, with Fit 5 corresponding to the whole sample and Fit 6 to the best-observed SNe. 
As is seen, these fits still prefer a value of Ry ~ 1.5. 

Bearing in mind the slight differences between photometric systems, the slope of the luminosity 
vs. decline-rate relation found for Fits 1-6 is in r easonable agreem ent with previously published 
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As shown in the left panels 
Phillips et al.l ()1999l ) is also consistent with the 



results (jHamuy et al.lll996al : lAltavilla et alJl2004l ; 
of Figure \T7\ the quadratic relationship found by 
data, although the CSP observations themselves are fit perfectly well by a simple linear relation. 

Next, we considered fits using the peak magnitudes in ugriVYJHK s , combining each band 
with a suitable color excess. The choice of this color excess is somewhat arbitrary. For the ugr and 
JHK S bands, color excesses involving the filter itself were used. However, when this is done for 
the i and Y filters, the resulting fits yield significantly higher dispersions than are obtained when 
using E(B — V). We attribute this to the poor quality of the template fits in the i and Y bands 
(see § I2.3p . which not only introduces an error in the peak magnitude, but also in the color excess 
if the latter involves either of these bands. 

Fits 7-15 in Table[9]give the results for all bands other than B using the sample of best-observed 
SNe. The best-fit values of Ry range from 1.1 to 1.8 for ugriYJ. For H and K s , somewhat larger 
values in the range 2.4-2.7 are obtained, but with such large uncertainties that they are compatible 
with both low and standard values of Ry. A weighted average of the results for Fits 7-15 yields 
Ry = 1.38±0.04, which is consistent with the results obtained for B vs. E(B — V) (see Fit 4). The 
scatter in the corrected magnitudes for Fits 7-15 is in the range 0.12-0.16 mag, with the intrinsic 
components of the dispersion amounting to as small as 0.04 mag and as large as 0.13 mag. 



4.3. SNe la as Standard Candles in the Near-Infrared 



The set of NIR light curves presented here constitutes a significant contribution to the avail- 
able data in the literature. The advantages of using NIR observations of SNe la to determine 
distances were demonstrated by iKrisciunas. Phillips. &: Suntzeff I ([2004) and, more recently, by 



At maximum light, the B — V colors of SNe la are close to zero, except for the most heavily reddened events. By 
definition, at B — V = 0.0 mag, magnitudes in the natural system are identical to those that have been color-corrected. 
Hence, comparing the CSP natural magnitudes at maximum light with those in systems where a color-term correction 
has been applied is a reasonable proposition. 



-22- 



Wood-Vasev et al.l (]2008l ). Apart from being less affected by dust extinction and scattering - 
which considerably reduces the effect of uncertainties in the color excesses and reddening law — the 
peak luminosities of SNe la in the NIR present a shallower dependence on decline rates than their 
optical counterparts. 

Indeed, previous studies have shown no clear evidence for a correlation be tween luminosity 



and d ecline rate in the JHK S bands except for the very fastest declining SNe (jKrisciunas et al 



2004bl ) . The improved quality and coverage of the CSP NIR light curves allow a weak dependence 
of the luminosity on the decline rate to be discerned in the J band. Fit 13 of Table [9] yields a slope 
of bj = 0.58 ± 0.09 over 0.7 < Amis(B) < 1.7 mag, and is shown in panel (c) of Figure [T71 In the 
H band, the evidence is also suggestive for a weak correlation between luminosity and decline rate, 
althou gh the significance o f the measured slope is only ~ 2a (see Fit 14 of Tabled]). In a future 
paper (jKattner et al.ll2009l ). we will use a larger sample of CSP SNe to test the reality of these 
relations. 

To compare with previous studies, we compute the dispersion in absolute peak magnitudes 
without any correction for decline rate. Table [10] shows the results of weighted averages of the 
absolute peak magnitudes in YJHK S . These averages correspond to the sample of best-observed 
SNe, correcting for extinction using the values of E(V-X) and R Y X Z listed in Fits 12-15 of Tabled 
For comparison, averages were also taken limiting the sample to those SNe with low/moderate 
reddening (i.e., excluding SN 2005A and SN 2006X), and assuming two different values of R\ Z '- 
zero, which is equivalent to ignoring host-galaxy extinction corrections, and the values of RX^ Z that 
correspond to Ry = 3.1. 

The most striking implication of Table [T0l is that, when the highly reddened SNe are excluded 
from the sample, the absolute magnitudes change only very slightly depending on the value of Ry. 
This illustrates the great advantage of working in the NIR, where lack of knowledge of the exact 
value of Ry has little influence, and where dust- extinction corrections can essentially be ignored for 
all but the most heavily reddened SNe. 

Table [TOl shows that the dispersions in the absolute magnitudes are 0.18-0.21 mag in the J 
and H bands, depending only slightly on the assumed value of Ry. These are to be compared 
with the rms dispersions of 0.12 mag and 0.16 mag, respectively, obtained when the absolute 
magnitud es are corrected f or the decline-rate dependence (cf. Fits 13 and 14 of Tabled]). For 
reference. iKrisciunas et al.1 (]2004bl ) obtained rms valu es (uncorrected for decline rate) of 0.13 mag 
in J and 0.15 mag in H for a sample of 22 SNe la, and IWood-Vasey et al.l (|2008l ) found dispersions 
of 0.33 mag in J and 0.15 mag in H for 21 SNe. 

In Table [11] the average absolute magnitude s at maximum of the be st-observed subsample 
of CSP SNe are compared with those obtained by IKrisciunas et al.l (|2004bl ) . The values in J and 
H are consistent at the ~ 2a level. In this same table, a comparison is als o attem pted with the 
absolute magnitudes at the epoch of B maximum given by I Wood- Vasey et al.1 (|2008l ). For J and H, 
these were derived using all the measurements available for the CSP SNe inside a bin of [—2.5,2.5] 
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rest-frame days with respect to B maximum, converting these to absolute magnitudes as explained 
below. A total of 33 and 28 data points entered into the calculations for J and H , respectively. For 
the K s band, the relative scarcity of observations forces the use of a slightly larger bin size of [—4,4] 
days in order to include 14 data points in the average. Again, the results in J and H are consistent 



with those of lWood-Vasev et al.l (|2008l ) at only the ~1-2ct level. In both cases, we suspect that the 
relatively poor agreement is due more to differences in the methods used to fit and combine the 
data rather than differences in the data themselves. 

Figure [18] shows the absolute iYJH light curves for the sample of best-observed CSP SNe. 
The time axis is plotted as rest-frame days since maximum light in B. The Hubble- flow distances 
of equation ([5]) along with the Cepheid and SBF distances given in Table [7] are used to convert the 
observed magnitudes to absolute values. K-corrections (see § [2]), Galactic-extinction corrections 
with Ry al = 3.1, and host-galaxy extinction corrections using the measurements of color excesses 
E{V — X) and the Rv Z parameters obtained from Fits 11-14 for iYJH have also been applied. 
For the fast-declining SNe 2005ke and 2006mr, zero host-galaxy extinction is assumed. 

Ignoring for the moment the two fast-declining SNe 2005ke and 2006mr, Figure [18] indicates 
that the smallest dispersion in absolute magnitude occurs around the time of the primary maximum 
in each filter. By the time of the minimum between the primary and secondary maximum, the 



dispersion increases in all bands. This observation is at odds with the suggestion by iKasenl ((2006) 
that the minimum of the J-band light curves of SNe la should show a low dispersion in luminosity. 
At epochs later than the second maximum, the data in iYJH show increasing spread due to 
differences in the decline rates among the SNe (see also § [2] and Figure [2]) . 

The two fast-declining events (SN 2005ke and SN 2006mr) do not reach maximum brightness 
until ~5 days after the epoch of B maximum. Interestingly, Figure [18] shows that at this epoch, 
their J-band luminosities are similar to those of the rest of the SNe which are already declining 
from the first maximum. 



4.4. Combined Hubble Diagram 

The fits presented in § 14.21 can be used to place the SNe with 0.7 < Amis(-B) < 1.7 mag 
in a Hubble diagram. Distance moduli were computed for each SN in ugriBVYJH using the 
parameters given in Table [9] for the fits to the best-observed SNe (Fits 4 and Fits 8-14). Note that 
Fit 7 involving V and E(B — V) was not included because it is mathematically equivalent to that 
involving B and E{B — V) (Fit 4). The distance moduli for all available bands were then averaged 
for each SN. The full correlation matrix among distance moduli from all bands was estimated in 
order to weight the average. Correlation might be expected to arise from the fact that all fits 
involve the same decline-rate parameter Ami^(B), and in several cases, the same bands in the 
color-excess term. These averaged values are plotted versus the luminosity distance, di, in the 
top half of Figure [19] A combined scatter of 0.11 mag, or ~5% in distance, is found for the 23 
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best-observed SNe that are included in the diagram. 

In the bottom panel where the residuals from the fit are plotted, it may be seen that the 
dispersion decreases with distance, implying that it is due predominantly to the peculiar velocities 
of the individual host galaxies. Plotted for reference as d otted lines is the spread predicted by the 



la dispersion of 382 km s , as derived by I Wang et al.l ( 2006 ) from an independent sample of 56 



SNe la. In general, the distance modulus residuals are consistent with the latter value, although 
it is interesting that the four most-distant SNe all have negative residuals (meaning that they are 
overluminous) . The most discrepant of these objects is SN 2004gu which, as discussed in C09, 
appears to have been spectroscopically similar to the peculiar SN la 2006gz, a slow-declining, over - 



luminous SN la that may have resulted from the merger of two white dwarfs ([Hicken et al.l 120071 ) . 
Both SNe displayed nearly identical decline rates (Amis(S) ~ 0.7), and very similar luminosities 
(Mb ~ —19.6 mag). The other three SNe (2005ag, 2005ir, and 2006py) are all slow-decliners 
(0.86 < Ami^,(B) < 1.02 mag) and, therefore, luminous events. In principle, the luminosity vs. 
decline-rate relationship should correct for this. Nevertheless, SN 2005ag was discovered by the 
LOSS, and must have been at the limit of detection. Hence, Malmquist bias could explain its over- 
luminous nature. However, SNe 2005ir and 2006py were discovered by the much deeper SDSS II 
survey which must have negligible bias at these redshifts. SN 2004gu had a moderately large color 
excess of E(B — V) max = 0.20, but the other three events had color excesses very close to the 
median value of E(B — V) max = 0.05 for the full subset of best-observed SNe. Hence, it is unlikely 
that the negative residuals in Figure [19] can be ascribed to inaccurate reddening corrections. Since 
the points for SNe 2005ag, 2005ir, and 2006py all lie within la of the expected velocity spread due 
to peculiar velocites, we consider small-number statistics to be the most likely explanation. Clearly 
it will be interesting to re-examine the Hubble diagram residuals once the full CSP sample has been 
produced. 

Interestingly, although the Hubble diagram in Figure [19] was produced by averaging the dis- 
tance moduli derived in each filter, the resulting dispersion is not much better than the dispersions 
obtained in the individual filters (cf. Table [9]). The explanation for this is found in Figure [20] where 
the distance-modulus residuals in ugriYJHK are plotted versus the residual in B. A significant 
correlation is observed between the distance- modulus residuals in one band vs. those in another, 
particularly in the ugri bands. For this reason, the combined residuals do not reduce significantly 
the scatter in the diagram as compared with the fits to the individual bands. 

We measured correlation coefficients, r, among residuals in Bugri to be 0.8-0.9, and 0.3- 
0.6 between any band and YJHK. Most likely the lower correlation in the NIR is attributable 
to the difficulties of measuring precise peak magnitudes in these bands as discussed earlier. We 
compared the observed correlation coefficients with those estimated for each SN based on the 
measurement uncertainties and the form of the model in equation ([7]). We found the expected 
correlation coefficients to be r < 0.05 for all cases, except between B and i bands where we 
estimate r ~ 0.6, and between B and Y where we expect r ~ 0.3. For these two cases, the high 
correlation may in part be due to the use of E(B — V) in the color-excess term. For all the rest of 
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the cases, the contribution to the correlations due to the fit model is negligible. We thus interpret 
the dispersion in Figure [19] to be mostly due to peculiar velocities. If this is correct, then the true 
precision of the SN distances is the dispersion about the relations in Figure [20"1 which for the ugri 
bands amounts to 0.06-0.09 mag in distance modulus, or 3-4% in distance. 



5. DISCUSSION AND CONCLUSIONS 

This paper presents a first analysis of the light curves of 34 SNe la followed by the CSP and 
released by C09. The high photometric precision and dense time sampling of the observations, 
especially in the optical bands, has allowed an in-depth examination of the general properties of 
the light curves of SNe la in the ugriBVY J H bands. Subsamples of well-observed SNe — 26 
in the optical and 9 in the NIR — were used to build a family of template light curves using a 
technique that allows one to interpolate among the template data, taking into account variations in 
the sampling, in a 3-D space parameterized by the epoch relative to maximum light, the magnitude 
relative to maximum, and the decline rate. 

The availability of observations covering a wide range of wavelengths (from u to K s ) offers 
the opportunity to make further progress on the difficult issue of correcting SN luminosities for 
extinction outside the Galaxy. However, depending on the approach taken, two quite different 
results are found. In the first case, a subsample of 15 SNe assumed to have suffered little or no 
extinction in their host galaxies was used to measure color excesses for the whole sample. By plotting 
the ratios of color excesses involving optical and NIR filters and comparing the values predicted 
by the CCM+O extinction law for interstellar dust, a value of the total-to-selective absorption 
coefficient of Ry ~ 1.7 was derived, which is significantly lower than the Galactic average of 3.1. 
However, this value is largely influenced by two very red objects in the sample, SNe 2005A and 
2006X. When the same calculations are repeated after excluding these two objects, a value of 
Ry = 3.2 ± 0.4 is found, in agreement with the Galactic average. 

An alternative way of estimating Ry is to express the absolute magnitudes of the SNe as a 
two-parameter function of the decline rate and color, and then to use the Hubble- flow distances (or 
Cepheid and SBF distances) of the S Ne to d e rive t he be st fit parameters throug h \ 2 minimization. 



This method was first employed by I Tripp! (119981 ) and iTripp Branch! (|1999l ). who related the 



absolute magnitude in B to the B — V color and derived values of Ry ~ 1—1.5. When this same 
technique is applied to the CSP sample of SNe la using the peak magnitudes in B and J bands 
and the B — V and V — J colors, a value of Ry ~l-2 is obtained, regardless of whether the heavily 
reddened events, SNe 2005 A and 2006X, are included or excluded. An equivalent analysis can also 
be carried out relating absolute magnitude to the decline rate and color excess. Performing such fits 
for all available bands {ugriBVY J H K s ) again leads to systematically low values of the CCM+O- 
law parameter (Ry 1.5) regardless of whether the two highly reddened SNe are included in the 
sample. 
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These conflicting results on the value of Ry reflect the variety of results found in the literature 
on this subject. Most early attempts to derive the average propertie s of the host-galaxy reddening 
indicated that Ry was unusually low (see iBranch &: Tammannlll992i . and references therein) . How- 
ever, the available light curves at the time were largely photographic and there was considerable 
additional uncertainty regarding the nature of the intri nsic colors of SNela at maximum light. The 
first modern study based on CCD data was made by iRiess et al.l (119961 ). who used the ratios of 
color excesses measured in B — V, V — R, and V — I for a sample of 20 SNe la to derive a value of 
Ry = 2.55 ± 0.30. This analysis differed from previous studies in that the intrinsic color variation 
as a function of decline rate was taken fully into acc o unt in deriving the c o lor excesses. Simila r 
valu es of Ry have been found bv lPhillips et al.l dl999h. lAltavilla et al.l |2004h . iReindl et al.l (|200a), 
and IWang et al.l (|2006l ). However, beginning with iTrippI (|1998l ). studies based on the technique of 
minimizing the dispersion in the Hubble diagram of SNe la have nearly uniformly led to va lues 



of Ry in the range of ~l-2 (e.g., iTripp &: Branch! Il999l ; lAstier et al. 



2006 



Conlev et al.ll2007l ). If 



there is a pattern to these results, including our own in this paper, it would seem that attempts to 
measure Ry via comparison of colors or color excesses tend to give larger values than the procedure 
of minimizing the scatter in the Hubble diagram with Ry tre ated as a free parameter. Nevertheless, 
there are exceptions to this rule such as the recent paper by lNobili Goobarl (120081 ). who studied 
the color evolution in U — B, B — V, V — R, and R — I of 69 SNe la with moderate reddening and 
obtained a value of Ry = 1.01 ± 0.25. 

Common sense suggests that at least some portion of the observed reddening of SNe la in 
spiral galaxies like the Milky Way must be due to dust in the interstellar medium of these galaxies. 
Although only a few studies have been carried out of the nature of the reddenin g law in external 
galaxies, these are generally consistent with Ry ga 3, albeit with a large dispersion (IGoudfrooij et al. 
1994 : IPatil et al.l I2QO7I : bstman. Goobar. & Mortselll I2OO8I ) . As shown in § l3XT| the ratios of the 
measured color excesses involving optical and NIR bands w i th res pect to E(B — V) are consistent 
with values of Ry ~ 3. As discussed by iKrisciunas et al .1 (|2_007) , inclusion of the NIR bands in 



this type of analysis improves considerably the accuracy with which Ry can be measured. Hence, 
it is tempting to conclude that we are measuring reddening produced by dust with very similar 
properties to that found in the interstellar medium of the disk of the Milky Way. It will be 
interesting to see if this finding holds up as many more SNe la with reliable optical— NIR color 
excesses are added to the CSP sample. 

If we are, indeed, observing host-galaxy reddening due to "normal" Galactic-type dust, then 
the fact that low values of Ry (~l-2) are obtained when it is treated as a free parameter in 
minimizing the dispersion in the Hubble diagram is puzzling. The implication is that there is an 
intrinsic dispersion in the colors of SNe la which is correlated with luminosity, but is independent 
of the decline rate. Our finding in § 13.4.11 that there is a small but systematic difference between 
color-excess measurements, E(B — V), made at late epochs using the Lira law and those derived 
from the maximum light magnitudes, and that this difference correlates with color and perhaps 
also with absolute magnitude, may well be related to this. 
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Very red SNe potentially allow the dust reddening law to be studied on an individual basis. 
Including the results for SN 2005 A given in this paper, there are now six SNe la for which a 
determination of Ry has been made from optical and NIR photometry. The results for these 
objects are summarized in Table [T2l A weighted mean of the six determinations gives Ry = 1.6, 
with a surprisingly small rms of 0.3. The relatively narrow range of decline rates is also noteworthy, 
averaging Ami5(B) = 1.2 w ith a dispersion of only 0.1. Ex cept for SN 2001el, the color excesses for 



these SNe are all very large. iJha. Riess. &: Kirshneri (|2007l ) found that the distribution of E(B — V) 



for SNe la is well approximated by an exponential function with r = 0.138 ± 0.023 mag. Hence, 
values of E{B — V) > 1 mag should be extremely rare, although it must be kept in mind that 
the sample of SNe analyzed by Jha et al. was culled from several sources with different selection 
criteria. Assuming that the Jha et al. distribution is correct, less than one SN for every thousand 
observed would be expected to have such a large amount of reddening — yet according to the 
Asiago Supernova Catalog^], the total number of nearby (z < 0.02) SNe la discovered from 1985 
through 2008 was only 358. The fact that these heavily reddened SNe are much more common 
than expected, combined with the remarkable similarity of decline rates and Ry values for the 
five objects in Table [12] with E(B — V) > 1 mag, leads us to speculate that these events actually 
represent a physically distinct subclass of SNe la. Hence, while extremely interesting in themselves, 
these objects like SNe 2005A and 2006X are not representative of the class of SNe la employed in 
cosmological studies. 



Wand (|2005l ) has suggested that the presence of circumstellar dust distributed in a shell around 
the SN can produce reddening compatible with the observed low values of Ry. This hypothesis is 
supported by our observations of SNe 2005A and 2006X. Specifically, we find that the red colors of 
these two SNe are better matched by a mo del where mult iple scattering of photons by circumstellar 



dust steepens the effective extinction law (jGoobarl 120081 ) . than by the standard CCM+O Galactic 
reddening law with a low value of Ry. 

Our results on the absolute peak magnitude calibrations for this new sample of CSP SNe 
confirm those obtained in previous studies, while extending the calibration to the opt ical uqri and 



the NIR YJHK S bands. Fits to the subsample of best-observed SNe employing the iTrippI (119981 ) 
model, which assumes the absolute magnitudes to be a two-parameter function of the decline rate 
and color, give rms dispersions of 0.12-0.15 mag. Remarkably, the same fits are found to apply 
equally well to either highly-reddened objects such as SNe 2005A and 2006X, or fast-declining, 
intrinsically-red events like 2005ke and 2006mr. The simplicity of this method is appealing. We 
have also carried out fits which assume the absolute peak magnitudes to be a function of the 
decline rate and the measured color excesses. These give quite similar results to those obtained 
with the Tripp model. Fits to a large number of filters and colors give dispersions in the range of 
0.12-0.20 mag for SNe with 0.7 < Ami 5 (B) < 1.7 mag. 

The quality of our NIR light curves allows, for the first time, detection of a weak dependence 
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of the J-band luminosity on decline rate. Combining the calibrations from all bands, we created 
a single Hubble diagram for the 23 best-observed SNe. Although the resulting scatter of 0.12 mag 
(6% in distance) is excellent, it does not represent a significant improvement over the precision 
obtained with individual filters. This is attributed to the fact that, for any particular SN, there 
is a significant correlation in the error in the distance moduli obtained using different filters. The 
source of these correlated errors appears to be the peculiar velocities of the SN host galaxies. If 
true, this implies that the derived distances to the SNe are actually precise to 3-4%, making SNe la 
competitive with even Cepheid variables as extragalactic distance indicators. 

Finally, the set of NIR light curves presented here confirms the advantages of working in this 
wavelength region for cosmological studies. Fit 8 of Table [8] and Fit 13 of Table [9] show that using 
either the V — J color or color excess to correct the J absolute magnitudes of the best-observed 
SNe (excluding the fast decliners) yields a dispersion of only 0.12 mag. Moreover, this result is 
insensitive to the exact form of the reddening law since the extinction corrections in the NIR are 
small. Note that the values derived for a r mSN of 0.02-0.04 mag suggest that the intrinsic dispersion 
in the J band may be quite small. This is obviously a preliminary result which must be treated 
with caution — however, we will soon have a sample of ~80 SNe la to test this. If confirmed, 
it implies that if the observational uncertainties can be reduced, there is much to be gained by 
extending the rest-wavelength coverage of future dark-energy experiments to include the J band. 
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A. Chi-Squared Fitting 

The problem of calibrating the low-redshift CSP SNe in a particular band X amounts to fitting 
a surface in a 3-dimensional parameter space defined by (1) a decline-rate parameter, Ami5(B); 
(2) a color parameter c, which could represent a pseudocolor at maximum, a color excess, or some 
other color indicator; and (3) the luminosity, Mx- For our analysis, we choose to fit this surface 
with a simple plane: 

Mx(Am 15 (B),c) = Mx{0) + b x [Am 15 (B)-l.l] + p x c. 

The Mx values themselves are derived from the observed magnitudes mx by a distance modulus 
/i which, for the closest objects, can be derived from independent data (e.g., Cepheids, SBF). For 
the more distant objects, we use Hubble's law, assuming Hq = 72 km s _1 Mpc -1 . Our model for 
the observed magnitudes is therefore (cf. equations ([6]) and (jTJ)) 



m x = M x (0) + 6 x [Ami 5 (fl)-l.l] + p x c + M . (Al) 

In the case where c is interpreted as a color excess, we use the intrinsic color relation derived 
in § 13.31 to compute 

c = E{Y - Z) = (my - m z ) - A - B [Am l5 (B) - 1.1]. (A2) 

Substituting this in equation (|Aip gives 

m x = (M x (0) - I3 X A ) + (b x ~ PxB ) [Am 15 (B) - 1.1] + (3 X c + fx. (A3) 

In the case where c is implemented as a pseudocolor (§ I4.ip . the same model can be used, simply 
setting Aq and Bq to zero. In this way, we see that the two treatments of c are mathematically 
equivalent. We also see very clearly that any uncertainties in the values of Aq and Bq are to be 
treated as systematic errors and should not be included in the statistical error budget for the color 
excesses. Rather, they should be included in the systematic error budget for the Mx(0) and bx- 

Given a set of N SNe with measured values {/ij; mxi', Am\§{B)i; q} for i = 1,...,N, we 
construct the x 2 function 

N , _ , 2 

x\M x {%bx,fix,o SN ) = * 2 . (A4) 

,- = i a i + °SN 



where af are the measurement variances, and asN is an additional constant term to account for 
any possible intrinsic disp ersion of the SN data which has not been taken into account by the model 
(see iTremaine et al.ll2002l ). The distance moduli, /ij, are computed using equation ([5]). 
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A.l. Variances 

In computing the denominator of the x 2 function, care must be taken in estimating not only 
the variances of the data, but also the covariances between the different observables. We outline 
below how these variances and covariances are estimated. For each supernova, the total variance 
is (we omit the subscript i for clarity) 

a 2 = a\fi) + a 2 (m x ) + b 2 x a 2 (Am l5 (B)) + (3 2 x a 2 (c) 
- 2 6xcov[m x ,Ami 5 (B)] - 2 (5 X cov [m x , c] 

+ 2b x p x cov[Ami 5 (B),c], (A5) 

where <r 2 (q) denotes the variance in quantity q, cov[p, q] denotes the covariance between quantity 
p and q, and we have assumed cov[//, q] = 0. Furthermore, we assume the peak magnitudes in 
two different bands are uncorrelated (cov[m x , ttiy] = 5 X y(y 2 {m x )). The three covariance terms in 
equation (|A5[> are computed as follows. 



For cov[mx, Amis (5)], we recall that Ami^{B) = rusty = 15) — m B ty = 0) and so 

cov[mx, Ami5(-B)] = cov[m x , tub (t = 15)] — cov[m x , rusty = 0)] 

= SxBcr 2 (m x ), (A6) 

where we have assumed cov[m x , rusty = 15)] = 0. Note that we have neglected any possible 
correlated er ror introduced by the subtraction of the host-galaxy image in the B band. As 



discussed by IContreras et al.l (120091 ) , such errors are less than the quoted statistical errors of 



the photometry. 

• For cov [m x , c] we have, from equation (|A2|) , 

cov[mx,c] = cov[m x , (my — mz)] — Bo cov[mx, Amis(B)] 

= ($xy - Sxz + B 5 XB )(T 2 (m x ). (A7) 

• For cov[Ami5(B), c], we also have from equation (IA2I) 

°"(Am ls (B),c) = cov[Am 15 (B),(m Y - m z )] - cov[Am 15 (B),B Ami 5 (B)] 

= {5 Z B-SYB)a 2 m B - B a 2 (Am 15 {B)). (A8) 



Substituting equations (jA6j), (|A7j) . and (|A8j) into equation (|A5[) for a 2 we obtain the final 
expression for the measurement variances, 

a 2 = a 2 (/x) + [1 + 25 XB b x - 2(5 X {5 XY - S xz + B 5 XB )] a 2 (m x ) + 

+ [b\ - 2b x (3 x B ] a 2 (Am 15 (B)) + f3 2 x a 2 + 2b x (3 x {5 ZB - 5 YB )<J 2 (m B ), (A9) 

where we recall that Bq = for the case when c represents a pseudocolor. 
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A.2. Solution 

The y 2 minimizati on is performed by analytically marginalizing equation (|A4|) over Mx(0) 



(see iGoliath et al.l l200ll ) at each point in a grid of values of {bx', Px}, followed by a search for the 
minimum value of x 2 ■ This process is repeated for various values of cjsn — usually between and 
0.3 mag — in order to find the value <jsn that produces a minimum reduced x 2 of Xmin/ — 3) = 1. 

Once ctsn is found, it is used for computing the probability density p(Xmarg) a ^ each point in 
the grid of {bx; fix} values. This probability is given by 



Finally, the expectation value of each parameter, 9, in the grid of 9 = {Mx(0); bx', fix} is computed 
as the moment 

* _ fo P (e)d9 



Jp(9)d9 

and the variances in the fit parameters, as, are estimated as the second moment 



(All) 



_ J(9-9) 2 p(9)d9 



a i = ^ — . (A12) 

fp(9)d9 



B. Deriving Ry 



Here we explain how we compute the coefficients ax and bx for each band ugriBVY JHK S 
that allow us to derive Ry from the slopes E(V — X) max / E(B — V) ma _ x of § 13.4.11 and from the fit 
parameters R x z of § I4.2L 



The reddening law is assumed to be a function of wavelength of the form given bv lCardelli. Clayton, fc Mathis 



diasa): 



A x b x 



(Bl) 



We simulate different amounts of extinction Ay on the typical spectrum of a SN la at maximum 
light by multiplying equation (IBID into t he spectrum . The spectrum we use is the template SN la 
spectrum of epoch 20 from iHsiao et al.1 (|2007l . 120091 ) which corresponds to maximum light in B. 
We perform the multiplication for various values of Ry, computing synthetic photometry from the 
resulting spectra using the transmission functions of the CSP filters described in C09. The ob- 
served absorptions, Ax and Ay, are then given by the difference between the synthetic magnitudes 
obtained from the reddened template spectrum and those from the original unreddened one, for 
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bands X and V. We finally fit a linear relation between Ax and Ry as in equation ([BT]), and thus 
derive ax and bx for all bands ugriBVY JHK S . 



We have performed these calculations for the original ICardelli. Clayton. Mathisl (|l989l ) 
(CCM) reddening law, and also with the modifications introduced by lO'Donnelll (| 19941 ) (CCM+O). 
Columns 2-5 of Table [13] list the values of ax and bx obtained for each case. As will be explained 
below, these coefficients serve to convert an observed reddening measurement, such as those in 
§ 13.4.11 and § 14.21 to a value of Ry. We note that this would be a true value of Ry — that is, the 
parameter defining the reddening law of equation (|Bip . 



If we instead consider the observed values of Ry derived from synthetic photometry as Ry° B = 
Ay** / E(B — V) obs , we obtain different coefficients ax and bx from the fits of Ax versus Ry> s . These 
coefficients are listed in columns 6 and 7 of Table [TBI for the case of the CCM+O law. 

The form ulae to convert observe d slopes E(V — X) masi / E(B — V) max to Ry are given in 
Appendix A of lKrisciunas et al.l (|2006l ) [see their equations (A5) and (A6)]. We briefly explain here 
how to derive Ry from fit parameters R x z of § I4.2L 

; From the definition of R x = A X /E(B - V), R x z can be expressed as 



r>YZ 
K x 



Rx 



Ry — Rz 



(B2) 



Substituting equation (|B1[) into this expression, and recalling that A\jAy is equivalent to 
R\/Ry, we derive 



Ry 



R Y x Z (by - b z ) - b x 



R\ Z (ay 



az) 



(B3) 



By error propagation of the uncertainty in R\ z ■, the uncertainty in Ry becomes 



(b Y - b z ) [R Y X Z (or - az) - ax] + (ay - az) [bx - R Y X Z 



{by - b z )\ 



[R\ z (a Y - a z ) - axY 



a nYZ. 



(B4) 
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Table 1. 









E(B - V) 


Amie(B) 


(m) 


Best 


SN 


•^hclio 


ZCMB 


Gal (mag) 


(mag) 


(mag) 


Obs.? 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


2004ef 


0.0310 


0.0298 


0.056 


1.36(01) 


35.46(03) 


yes 


2004eo 


0.0157 


0.0147 


0.108 


1.37(03) 


33.92(05) 


yes 


2004ey 


0.0158 


0.0146 


0.139 


0.98(02) 


34.06(05) 


yes 


2004gc 


0.0321 


0.0307 a 


0.206 


1.09(06) 


35.34(04) 


no 


2004gs 


0.0267 


0.0275 


0.031 


1.54(01) 


35.51(03) 


yes 


2004gu 


0.0459 


0.0469 


0.026 


0.80(04) 


36.41(10) 


yes 


2005A 


0.0191 


0.0184 


0.030 


1.17(03) 


34.46(06) 


yes 


2005M 


0.0220 


0.0230 


0.031 


0.85(01) 


35.13(03) 


yes 


2005W 


0.0089 


0.0080 


0.072 


1.16(03) 




yes 


2005ag 


0.0794 


0.0801 


0.041 


0.86(01) 


37.62(08) 


yes 


2005al 


0.0124 


0.0133 


0.055 


1.22(04) 


34.01(06) 


yes 


2005am 


0.0079 


0.0090 


0.054 


1.48(02) 




yes 


2005be 


0.0350 


0.0356 


0.034 


1.47(08) 


35.85(05) 


no 


2005bg 


0.0231 


0.0242 


0.032 


1.02(08) 


35.08(06) 


no 


2005bl 


0.0241 


0.0251 


0.029 


1.80(04) 




yes 


2005bo 


0.0139 


0.0150 


0.046 


1.30(08) 


33.85(06) 


yes 


2005el 


0.0149 


0.0149 


0.114 


1.29(01) 


33.92(05) 


yes 


2005eq 


0.0290 


0.0284 


0.074 


0.72(02) 


35.44(03) 


yes 


2005hc 


0.0459 


0.0450 


0.029 


0.90(01) 


36.51(05) 


yes 


2005iq 


0.0340 


0.0330 


0.022 


1.24(01) 


35.80(04) 


yes 


2005ir 


0.0764 


0.0753 


0.030 


0.93(11) 


37.59(06) 


yes 


2005kc 


0.0151 


0.0139 


0.132 


1.14(01) 


33.94(05) 


yes 


2005ke 


0.0049 


0.0045 


0.027 


1.76(01) 




yes 


2005ki 


0.0192 


0.0204 


0.032 


1.27(02) 


34.57(04) 


yes 


20051u 


0.0320 


0.0313 


0.026 


0.88(03) 


35.98(04) 


no 


2005na 


0.0263 


0.0268 


0.077 


0.95(02) 


35.24(05) 


yes 


2006D 


0.0085 


0.0097 


0.046 


1.42(01) 




yes 


2006X 


0.0052 


0.0063 


0.026 


1.09(01) 


31.05(05) 


yes 


2006ax 


0.0167 


0.0180 


0.050 


0.99(02) 


34.30(04) 


yes 


2006bh 


0.0109 


0.0105 


0.026 


1.41(01) 


33.25(07) 


yes 


2006eq 


0.0495 


0.0484 


0.048 


1.66(06) 


36.54(02) 


no 


2006gt 


0.0448 


0.0437 


0.037 


1.66(03) 


36.49(02) 


yes 


2006mr 


0.0059 


0.0055 


0.021 


1.82(02) 




yes 


2006py 


0.0579 


0.0567 


0.065 


1.02(03) 


36.88(02) 


yes 
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a Value of zcmb corresponding to the systemic recession velocity of the 
group ARP 327 

Note. - Columns: (1) SN name; (2) Heliocentric redshift from the 
NASA/IPAC Extragalactic Database (NED); (3) Redshift in the frame 
of the 3K CMB (NED); (4) Galactic reddening (NED); (5) Observed 
Ami5(-B); (6) Average distance modulus from fits of § 14.21 (7) Does the SN 
belong to the "best-observed" sample?. 
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Table 2. Light-curve parameters 







^max 






m max 




Ami 


Epoch 




Fit 


Band 


-2, 


453, 000 




(mag) 




(mag) 


Range 


Npts. 


Method 


(1) 










(3) 




(4) 


(5) 


(6) 


(7) 
















SN 2004ef 










B 


264 


44+0 


05 


17 


059±0 


005 


0.721±0.014 


-8.5 


+63.2 


44 


Spline 


V 


266 


07±0 


03 


16 


900+0 


002 


1.120±0.008 


— 10.1 


+61.6 


48 


Spline 


u 


263 


01±0 


14 


17 


791±0 


012 


1.596±0.068 


—6.2 


+56.8 


41 


Spline 


9 


264 


82±0 


06 


16 


964±0 


003 


1.077±0.009 


—8.9 


+62.8 


50 


Spline 


r 


266 


34±0 


07 


16 


862±0 


003 


1.625±0.009 


— 10.3 


+61.3 


47 


Spline 


i 


261 


21±0 


08 


17 


415±0 


005 


1.199±0.013 


—5.3 


+66.3 


47 


Spline 


Y 


261 


25±0 


03 


17 


209±0 


034 




+8.8 


+69.5 


7 


Tempi. 


J 


261 


63±0 


03 


17 


230±0 


080 




+8.3 


+20.6 


5 


Tempi. 


K 


262 


47±0 


03 


17 


470±0 


473 


SN 2004eo 


+7.5 


+19.8 


6 


Tempi. 


B 


278 


78±0 


20 


15 


500±0 


008 


0.725±0.019 


— 12.0 


+61.8 


39 


Spline 


V 


280 


45±0 


11 


15 


358±0 


013 


1.070±0.021 


-13.7 


+60.1 


39 


Spline 


u 


277 


37±0 


12 


16 


064±0 


015 


1.786±0.022 


— 10.7 


+43.5 


37 


Spline 


9 


279 


23±0 


17 


15 


378±0 


014 


1.107±0.021 


— 12.5 


+61.4 


38 


Spline 


r 


280 


93±0 


27 


15 


358±0 


009 


1.522±0.019 


— 14.2 


+59.7 


39 


Spline 


i 


275 


88±0 


17 


15 


781±0 


010 


1.029±0.020 


—9.2 


+64.7 


39 


Spline 


Y 


274 


78±0 


07 


15 


923±0 


007 




—5.1 


+31.3 


15 


Spline 


J 


274 


67±0 


02 


15 


537±0 


006 




-5.0 


+62.9 


17 


Spline 


H 


274 


71±0 


17 


15 


791±0 


008 




— o.U 




18 


Spline 


K 


275 


66±0 


20 


15 


592±0 


010 


SN 2004ey 


—5.9 


+30.4 


14 


Spline 


B 


304 


39±0 


10 


15 


260±0 


005 


0.466±0.007 


—8.7 


+47.4 


30 


Spline 


V 


305 


70±0 


09 


15 


215±0 


006 


0.892+0.011 


— 10.0 


+46.1 


32 


Spline 


u 


303 


24±0 


08 


15 


675±0 


Oil 


1.277±0.015 


— 1.0 


+48.5 


25 


Spline 


9 


304 


73±0 


11 


15 


214±0 


003 


0.748±0.005 


-9.1 


+47.1 


26 


Spline 


r 


305 


95±0 


13 


15 


236±0 


004 


1.149±0.007 


-10.3 


+45.9 


32 


Spline 


i 


301 


18±0 


07 


15 


761±0 


003 


0.680±0.008 


-5.6 


+50.6 


32 


Spline 


Y 


299 


81±0 


04 


15 


785±0 


018 




+1.8 


+57.7 


12 


Tempi. 


J 


301 


06±0 


04 


15 


504±0 


050 




-0.5 


+56.4 


15 


Tempi. 


H 


300 


63±0 


04 


15 


743±0 


022 




+1.0 


+58.9 


16 


Tempi. 


K 


302 


05±0 


04 


15 


589±0 


046 


SN 2004gc 


-0.4 


+55.4 


14 


Tempi. 


B 


324 


57±0 


26 


17 


467±0 


028 




+5.4 


+75.4 


27 


Tempi. 


V 


325 


85±0 


26 


17 


136±0 


024 




+4.1 


+74.1 


28 


Tempi. 


9 


325 


04±0 


26 


17 


328±0 


026 




+4.9 


+74.9 


26 


Tempi. 


r 


325 


78±0 


26 


17 


039±0 


025 




+4.2 


+74.1 


28 


Tempi. 


i 


321 


55±0 


26 


17 


446±0 


026 


SN 2004gs 


+8.5 


+78.5 


26 


Tempi. 


B 


356 


18±0 


04 


17 


260±0 


005 


0.854±0.022 


-3.2 


+96.8 


51 


Spline 


V 


358 


64±0 


18 


17 


072±0 


005 


1.368±0.025 


-5.7 


+100.2 


55 


Spline 


u 


354 


98±0 


12 


17 


983±0 


018 


1.756±0.087 


-2.1 


+35.7 


23 


Spline 


9 


356 


79±0 


04 


17 


138±0 


004 


1.294±0.034 


-3.8 


+102.1 


51 


Spline 


r 


358 


19±0 


21 


17 


046±0 


005 


1.733±0.018 


-5.2 


+100.7 


52 


Spline 
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m™ 


Am t 


Epoch 




Fit 


Band 


-2,453,000 


(mag) 


(mag) 


Range 


Npts. 


Method 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 



i 


355 


76±0 


05 


17 


449±0 


007 


1.469±0.214 


-3.0 


+108.7 


50 


Tempi. 


J 


354 


61±0 


05 


17 


323±0 


017 


SN 2004gu 


+8.4 


+60.6 


19 


Tempi. 


B 


361 


91±0 


17 


17 


554±0 


022 




-0.1 


+51.3 


25 


Tempi. 


V 


363 


26±0 


78 


17 


380±0 


009 




-5.2 


+45.5 


26 


Spline 


u 


360 


54±0 


17 


18 


096±0 


038 




+1.4 


+35.9 


14 


Tempi. 


9 


362 


90±0 


33 


17 


417±0 


008 




-4.8 


+45.8 


26 


Spline 


r 


364 


77+0 


21 


17 


340±0 


005 




-6.6 


+44.0 


29 


Spline 


i 


359 


83±0 


17 


17 


916±0 


013 


SN 2005A 


-2.1 


+53.5 


27 


Tempi. 


B 


380 


17±0 


12 


18 


262±0 


008 


0.560±0.030 


-3.5 


+29.8 


23 


Spline 


V 


382 


72±0 


71 


17 


235±0 


007 


0.956±0.018 


-6.0 


+48.9 


32 


Spline 


u 


378 


85±0 


72 


19 


933±0 


043 




-2.2 


+11.5 


10 


Spline 


9 


381 


92±0 


35 


17 


830±0 


007 


0.847±0.026 


-5.2 


+49.6 


36 


Spline 


r 


381 


50±0 


17 


16 


844±0 


007 


1.283±0.015 


-4.7 


+50.1 


36 


Spline 


i 


376 


60±0 


29 


17 


184±0 


014 




+0.1 


+57.0 


35 


Tempi. 


Y 


374 


36±0 


29 


16 


398±0 


026 




+11.4 


+60.3 


17 


Tempi. 


J 


375 


62±0 


29 


16 


217±0 


024 




+10.1 


+56.0 


17 


Tempi. 


H 


375 


09±0 


29 


16 


404±0 


034 


SN 2005M 


+10.7 


+59.6 


16 


Tempi. 


B 


405 


84±0 


07 


16 


046±0 


004 


0.393±0.010 


-7.9 


+73.1 


56 


Spline 


V 


407 


88±0 


16 


16 


011+0 


003 


0.838±0.014 


-9.9 


+71.1 


58 


Spline 


u 


404 


46±0 


15 


16 


359±0 


018 


1.158±0.026 


-7.5 


+52.0 


47 


Spline 


9 


406 


28±0 


08 


15 


962±0 


004 


0.662±0.015 


-8.4 


+72.6 


62 


Spline 


r 


408 


64±0 


14 


16 


102±0 


006 


0.994±0.011 


-10.7 


+70.3 


59 


Spline 


i 


402 


98±0 


11 


16 


575±0 


004 


0.688±0.013 


-5.1 


+75.9 


59 


Spline 


Y 


402 


11±0 


14 


16 


610+0 


014 


0.071±0.021 


-4.2 


+58.3 


29 


Spline 


J 


402 


80±1 


44 


16 


495±0 


048 


1.650±0.054 


-4.9 


+57.6 


32 


Spline 


H 


401 


92±1 


93 


16 


717±0 


019 


0.216±0.024 


-4.1 


+59.3 


18 


Spline 


K 


403 


22±0 


02 


16 


403±0 


036 


SN 2005W 


-5.6 


+60.7 


15 


Tempi. 


B 


412 


08±0 


11 


14 


649±0 


005 


0.573±0.024 


-7.5 


+17.3 


20 


Spline 


V 


413 


10±0 


36 


14 


235±0 


009 




-8.5 


+16.3 


21 


Spline 


u 


410 


86±0 


16 


15 


014±0 


007 


1.277±0.119 


-6.3 


+16.5 


18 


Spline 


9 


412 


55±0 


14 


14 


330±0 


009 


0.791+0.027 


-7.9 


+16.8 


20 


Spline 


r 


413 


88±0 


53 


14 


196±0 


009 




-9.3 


+15.5 


21 


Spline 


i 


408 


44±0 


43 


14 


637±0 


009 


SN 2005ag 


-3.9 


+20.9 


19 


Spline 


B 


413 


53±0 


13 


18 


593±0 


006 




-0.8 


+73.2 


44 


Tempi. 


V 


415 


14±0 


13 


18 


550±0 


006 




-2.5 


+71.4 


44 


Tempi. 


9 


413 


98±0 


13 


18 


551±0 


005 




-1.3 


+72.6 


44 


Tempi. 


r 


415 


70±0 


13 


18 


631±0 


005 




-3.1 


+70.8 


46 


Tempi. 


i 


411 


23±0 


13 


18 


955±0 


020 




+1.7 


+75.6 


44 


Tempi. 


Y 


408 


79±0 


13 


19 


325±0 


021 




+10.9 


+45.7 


10 


Tempi. 
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m™ 


Am t 


Epoch 




Fit 


Band 


-2,453,000 


(mag) 


(mag) 


Range 


Npts. 


Method 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 



J 


410.18±0.13 


19.255±0.037 




+9.4 


+44.2 


10 


Tempi. 


H 


409.65±0.13 


19.067±0.055 


SN 2005al 


+10.0 


+44.8 


10 


Tempi. 


B 


430.84±0.30 


15.084±0.006 


0.630±0.009 


-1.1 


+80.9 


35 


Spline 


V 


431.57±0.38 


15.076±0.005 


1.066±0.010 


-1.8 


+80.2 


36 


Spline 


u 


429.12±0.11 


15.530±0.014 


1.565±0.017 


+0.6 


+84.7 


35 


Tempi. 


9 


430.65±0.15 


15.033±0.004 


0.935±0.009 


-0.9 


+81.1 


35 


Spline 


r 


431.91±0.48 


15.122±0.007 


1.573±0.012 


-2.1 


+79.9 


36 


Spline 


i 


427.98±0.11 


15.729±0.009 




+1.8 


+85.9 


35 


Tempi. 


Y 


426.43±0.11 


15.677±0.013 




+7.5 


+36.8 


12 


Tempi. 


J 


427.05±0.11 


15.531±0.015 




+6.9 


+36.2 


13 


Tempi. 


H 


426.67±0.11 


15.808±0.017 




+7.3 


+36.6 


13 


Tempi. 


K 


428.05±0.11 


15.720±0.064 


SN 2005am 


+5.9 


+35.2 


9 


Tempi. 


B 


436.86±0.07 


13.857±0.003 


0.775±0.007 


-4.3 


+75.0 


37 


Spline 


V 


438.45±0.12 


13.760±0.003 


1.256±0.007 


-5.8 


+73.4 


37 


Spline 


u 


435.63±0.05 


14.408±0.005 


1.750±0.014 


-3.0 


+76.2 


37 


Spline 


9 


437.28±0.03 


13.780±0.002 


1.152±0.008 


-4.7 


+74.6 


38 


Spline 


r 


438.43±0.10 


13.742±0.005 


1.741±0.009 


-5.8 


+73.4 


38 


Spline 


i 


435.36±0.08 


14.254±0.010 




-2.8 


+77.7 


36 


Tempi. 


J 


434.10±0.08 


14.009±0.014 




-0.4 


+28.7 


12 


Tempi. 


H 


434.26±0.08 


14.132±0.011 




-0.6 


+28.5 


14 


Tempi. 


K 


434.50±0.08 


14.135±0.043 


SN 2005be 


-0.8 


+28.2 


14 


Tempi. 


B 


460.59±0.28 


17.023±0.035 




+7.6 


+32.3 


13 


Tempi. 


V 


462.12±0.28 


16.974±0.030 




+6.0 


+30.7 


13 


Tempi. 


u 


459.00±0.28 


17.364±0.042 




+9.2 


+25.7 


9 


Tempi. 


9 


461.12±0.28 


16.950±0.032 




+7.0 


+31.7 


12 


Tempi. 


r 


461.90±0.28 


17.054±0.028 




+6.2 


+30.9 


12 


Tempi. 


i 


458.83±0.28 


17.630±0.020 


SN 2005bg 


+9.4 


+34.1 


12 


Tempi. 


B 


470.70±0.33 


16.010±0.042 




+1.0 


+42.9 


16 


Tempi. 


V 


471.89±0.33 


15.941±0.044 




-0.3 


+41.7 


16 


Tempi. 


u 


469.24±0.33 


16.400±0.037 




+3.5 


+22.9 


13 


Tempi. 


9 


471.13±0.33 


15.997±0.043 




+1.6 


+42.5 


15 


Tempi. 


r 


472.08±0.33 


16.082±0.046 




+0.6 


+41.5 


15 


Tempi. 


i 


467.74±0.33 


16.676±0.052 


SN 2005bl 


+5.1 


+45.9 


15 


Tempi. 


B 


484.25±0.36 


18.336±0.031 




-7.5 


+27.6 


12 


Spline 


V 


485.39±0.11 


17.753±0.017 




-8.6 


+26.5 


12 


Spline 


u 


481.42±0.67 


19.784±0.076 




-4.7 


+9.9 


10 


Spline 


9 


484.03±0.53 


18.129±0.057 




-7.3 


+27.9 


12 


Spline 


r 


486.31±0.12 


17.694±0.004 




-9.5 


+25.7 


12 


Spline 


i 


486.50±0.09 


17.791±0.006 




-9.6 


+25.5 


12 


Spline 
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m™ 


Am t 


Epoch 




Fit 


Band 


-2,453,000 


(mag) 


(mag) 


Range 


Npts. 


Method 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 



SN 2005bo 



B 


479 


76+0 


08 


15 


845±0 


009 




-1.0 


+32.6 


9 


Spline 


V 


480 


80±0 


14 


15 


515+0 


015 




-2.0 


+31.5 


10 


Spline 


u 


477 


62±0 


14 


16 


463±0 


036 




+1.1 


+14.3 


8 


Tempi. 


a 


480 


29±0 


22 


15 


701+0 


007 




-1.5 


+32.0 


9 


Spline 


r 


480 


95±0 


28 


15 


457±0 


010 




-2.2 


+31.4 


9 


Spline 


i 


476 


50±0 


14 


15 


790±0 


027 


SN 2005el 


+2.3 


+36.8 


9 


Tempi. 


B 


647 


24±0 


05 


15 


290±0 


008 


0.657±0.010 


-7.2 


+82.3 


25 


Spline 


V 


647 


89±0 


08 


15 


235±0 


005 


1.121+0.007 


-7.9 


+81.6 


24 


Spline 


u 


645 


52±0 


03 


15 


664±0 


006 


1.540±0.017 


-5.6 


+84.0 


24 


Spline 


9 


647 


34±0 


04 


15 


167+0 


010 


1.039±0.011 


-7.4 


+82.2 


24 


Spline 


r 


646 


89±0 


08 


15 


290±0 


009 


1.485±0.014 


-6.9 


+82.6 


24 


Spline 


i 


643 


75±0 


17 


15 


714+0 


006 


1.285±0.013 


-3.8 


+85.7 


25 


Spline 


Y 


642 


35±0 


03 


15 


514±0 


010 


0.377±0.018 


-1.5 


+72.3 


42 


Spline 


J 


642 


81±0 


13 


15 


584±0 


007 


2.443±0.044 


-1.9 


+71.8 


36 


Spline 


H 


642 


22±0 


13 


15 


642±0 


009 


0.995±0.032 


-1.3 


+72.4 


27 


Spline 


K 


644 


59±0 


07 


15 


572±0 


135 


SN 2005eq 


-2.9 


+37.6 


4 


Tempi. 


B 


654 


17+0 


14 


16 


591+0 


006 


0.321±0.011 


-3.3 


+94.7 


27 


Spline 


V 


655 


68±0 


25 


16 


453±0 


007 


0.793±0.017 


-4.8 


+93.3 


27 


Spline 


u 


652 


10+0 


09 


16 


914+0 


013 


0.894±0.344 


-1.3 


+62.4 


19 


Tempi. 


9 


654 


46±0 


16 


16 


496±0 


005 


0.562±0.008 


-3.6 


+94.5 


27 


Spline 


r 


656 


16±0 


13 


16 


485±0 


006 


0.966±0.013 


-5.2 


+92.8 


27 


Spline 


i 


652 


23±0 


09 


17 


032±0 


050 




-1.4 


+102.4 


27 


Tempi. 


Y 


649 


80±0 


09 


16 


902±0 


050 




-2.0 


+62.7 


28 


Tempi. 


J 


650 


64±0 


09 


16 


872±0 


029 




-2.9 


+61.8 


26 


Tempi. 


H 


650 


24±0 


09 


17 


090±0 


039 




-1.5 


+43.8 


11 


Tempi. 


K 


651 


28±0 


09 


16 


611±0 


073 


SN 2005hc 


+2.6 


+2.6 


2 


Tempi. 


B 


667 


14±0 


08 


17 


426±0 


004 


0.425±0.018 


-4.2 


+80.8 


24 


Spline 


V 


668 


88±0 


65 


17 


403±0 


006 


0.822±0.020 


-5.9 


+82.9 


25 


Spline 


u 


666 


28±0 


16 


17 


875±0 


009 


1.066±0.037 


-3.4 


+23.4 


11 


Spline 


9 


667 


69±0 


16 


17 


368±0 


005 


0.700±0.024 


-4.7 


+84.0 


23 


Spline 


r 


669 


01+0 


13 


17 


469±0 


005 


1.028±0.012 


-6.0 


+82.8 


23 


Spline 


i 


665 


13±0 


08 


18 


143±0 


008 




-2.5 


+79.1 


22 


Tempi. 


Y 


662 


73±0 


08 


17 


786±0 


050 




-1.1 


+50.2 


23 


Tempi. 


J 


664 


16±0 


08 


17 


936±0 


046 




-2.6 


+48.7 


22 


Tempi. 


H 


662 


92±0 


08 


17 


999±0 


103 


SN 2005iq 


-1.2 


+35.5 


13 


Tempi. 


B 


687 


78±0 


06 


16 


861±0 


007 


0.633±0.012 


-5.1 


+65.5 


20 


Spline 


V 


689 


70±0 


17 


16 


872±0 


010 


1.139+0.020 


-6.9 


+63.7 


21 


Spline 


u 


686 


49±0 


10 


17 


227±0 


010 


1.329±0.041 


-3.8 


+58.1 


17 


Spline 


9 


688 


16±0 


07 


16 


804±0 


007 


0.988±0.011 


-5.4 


+65.2 


20 


Spline 
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^max 






™max 
III 






Epoch 




Pit 


Band 


-2, 453, 000 




( ma g, 




(mag) 


Range 


Npts. 


A/T^tVinr] 
1V1C LIlLltl 






(2) 












(5) 


(°) 


(71 

(1) 


r 


689 


08±0 


17 


16 


957+0 


007 


1.482+0.018 


— O.O 


-(-04.0 


20 


Spline 


i 


684 


38±0 


09 


17 


421±0 


Oil 


1.111±0.025 


— l.o 


-j-Oo.o 


20 


Spline 


Y 


683 


67±0 


06 


17 


224±0 


050 




+2.0 


+50.7 


22 


Tempi. 


J 


684 


33±0 


06 


17 


402±0 


045 




I 1 Q 

+ 1.0 




19 


Tempi. 


H 


684 


21±0 


06 


17 


363±0 


050 


SN 2005ir 


-Li 7 




20 


Tempi. 


B 


685 


60±1 


26 


18 


518±0 


024 


0.483±0.030 


— Z.o 


/inn 
+4U.U 


12 


Spline 


V 


688 


25±1 


09 


18 


429±0 


026 


1.015±0.034 


— O.Z 




14 


Spline 


u 


683 


37±0 


23 


18 


889±0 


054 




i n ^ 
tU.o 


-^99 fl 


8 


Tempi. 


g 


685 


81±0 


61 


18 


497±0 


014 


0.776±0.024 


— o.U 


^9n q 


10 


Spline 


r 


687 


61±0 


33 


18 


512±0 


021 




— <±.o 


1 1Q 7 
TlO. 1 


10 


Spline 


i 


681 


96±0 


23 


19 


163±0 


053 


0.739±0.126 
SN 2005kc 


-un 7 

-+-U. i 




12 


Tempi. 


B 


697 


75±0 


06 


16 


033±0 


005 


0.560±0.012 


in 1 

— 1U. 1 


+ OU.O 


13 


Spline 


V 


699 


56±0 


12 


15 


742±0 


006 


1.028±0.012 


1 1 Q 




13 


Spline 


u 


696 


50±0 


16 


16 


704±0 


010 


1.465±0.034 


— O.O 


J-O'i 7 
T^O- 1 


11 


Spline 


g 


698 


39±0 


08 


15 


886±0 


005 


0.894±0.011 






13 


Spline 


r 


699 


92±0 


42 


15 


633±0 


005 




— 12 2 


^9fl ^ 


12 


Spline 


i 


695 


40±0 


05 


15 


924±0 


006 




7 £ 
— I .o 


7 

-rOA. 1 


13 


Spline 


Y 


694 


18±0 


10 


15 


649±0 


009 




— O.O 


J_9zl fl 


18 


Spline 


J 


694 


56±0 


11 


15 


494±0 


009 




— o.y 




16 


Spline 


H 


693 


99±0 


20 


15 


628±0 


039 


SN 2005ke 


— O.O 




16 


Spline 


B 


699 


16±0 


04 


14 


882±0 


005 


1.273±0.010 


Q A 

— y.4 




39 


Spline 


V 


701 


57±0 


07 


14 


185±0 


004 


1.657±0.011 


— 11.8 


+265.0 


39 


Spline 


u 


698 


05±0 


15 


16 


148±0 


008 


1.885±0.027 


-8.3 


+268.5 


36 


Spline 


g 


699 


92±0 


02 


14 


499±0 


003 


1.639±0.007 


-10.2 


+266.7 


39 


Spline 


r 


702 


20±0 


07 


14 


062±0 


003 


2.076±0.008 


-12.4 


+263.3 


40 


Spline 


i 


701 


80±0 


06 


14 


317+0 


004 


2.021±0.010 


-12.0 


+264.8 


39 


Spline 


Y 


701 


41±0 


17 


14 


076±0 


007 


1.390±0.018 


-9.6 


+87.7 


57 


Spline 


J 


700 


50±0 


06 


Id 


021±0 


007 


3.389±0.021 


-8.7 


+88.6 


54 


Spline 


H 


701 


76±0 


18 


1.3 


951±0 


007 


2.253±0.013 
SN 2005ki 


-9.9 


+87.4 


56 


Spline 


B 


705 


59±0 


11 


15 


705±0 


005 


0.641±0.009 


-9.5 


+152.9 


43 


Spline 


V 


706 


70±0 


09 


15 


672±0 


005 


1.143±0.010 


-10.6 


+151.8 


43 


Spline 


u 


702 


77±0 


10 


16 


050±0 


008 


1.336±0.032 


-6.8 


+80.4 


31 


Spline 


9 


706 


13±0 


07 


15 


629±0 


004 


1.063±0.009 


-10.1 


+152.4 


47 


Spline 


r 


707 


03±1 


21 


15 


737±0 


020 


1.651±0.022 


-11.0 


+151.5 


47 


Spline 


i 


701 


33±0 


23 


16 


208±0 


015 


1.225±0.019 


-5.4 


+157.1 


47 


Spline 


Y 


701 


75±0 


05 


16 


230±0 


013 


0.278±0.021 


-4.8 


+54.0 


22 


Spline 


J 


702 


08±1 


28 


16 


163±0 


072 


1.921±0.083 


-5.1 


+53.7 


21 


Spline 


H 


702 


03±0 


44 


16 


236±0 


027 


1.072±0.047 
SN 20051u 


-4.1 


+53.8 


17 


Spline 


B 


711 


90±0 


45 


17 


380±0 


018 




+9.0 


+90.5 


17 


Tempi. 



Table 2 — Continued 



t max m max Am( Epoch Fit 



Band 


-2,453,000 




(mag) 




(mag) 


Range 


Npts. 


Method 






(2) 






f3"l 




(A\ 


(5) 






v 


713 


27±0 


45 


17 


163+.0 


013 




J_7 f, 

-\- 1 .0 


Toy. i 


17 


Tempi. 




710 


36±0 


45 


17 


920±0 


058 




-h-LU.O 


-f-Zo.Z 


5 


Tempi. 




712 


35±0 


45 


17 


225±0 


014 




-ro.O 




10 


Tempi. 


r 


713 


73±0 


45 


17 


135±0 


Oil 






TOO.O 


16 


Tempi. 


i 


709 


23±0 


45 


17 


659±0 


015 


SN 2005na 


1 1 1 Q 
Tll.O 


i no o 

-t-yo.o 


16 


Tempi. 


B 


740 


02±0 


20 


16 


281±0 


009 


0.420±0.012 


-1.3 


+92.1 


26 


Spline 


V 


741 


70±0 


29 


16 


257±0 


006 


0.967±0.011 


o o 
— z.y 


i on a 


27 


Spline 


u 


738 


81±0 


17 


16 


442±0 


072 


1.369±0.231 


— U. 1 


J-A7 fi 


17 


Tempi. 


g 


740 


20±0 


31 


16 


235±0 


006 


0.740±0.008 


— 1.0 


T" (O.O 


26 


Spline 


r 


742 


25±0 


11 


16 


306±0 


006 


1.238±0.013 


— O.O 


1 on n 


27 


Spline 


i 


737 


51±0 


17 


16 


754±0 


050 




-1-1 9 




27 


Tempi. 


Y 


735 


52±0 


17 


16 


912±0 


017 




-\- -LO.U 


-U71 Q 


21 


Tempi. 


J 


736 


69±0 


17 


16 


738±0 


041 




J_1 7 A 


j_7n 7 


13 


Tempi. 


H 


736 


14±0 


17 


16 


975±0 


047 


SN 2006D 


Tlo.U 




16 


Tempi. 


13 


757 


95±0 


10 


14 


321±0 


006 


0.773±0.007 


(\ 1 
— 0. 1 


_|_1 1 -| Q 


39 


Spline 


V 


759 


24±0 


13 


14 


180±0 


004 


1.226±0.010 


-7.3 


+110.5 


41 


Spline 


u 


755 


54±0 


15 


14 


994±0 


013 


1.564±0.018 


— O. ( 


+ 04.4 


32 


Spline 


9 


758 


20±0 


08 


14 


225±0 


004 


1.097±0.006 


P. ^ 

— u.o 


4-1 11^ 
-(-±1 ±.o 


42 


Spline 


r 


758 


60±0 


09 


14 


185±0 


006 


1.670±0.010 


— U. ( 


+ 111 1 


42 


Spline 


i 


755 


23±0 


07 


14 


595±0 


005 


1.420±0.009 


— O.O 


-f-1 14.0 


42 


Spline 


Y 


754 


38±0 


05 


14 


650±0 


050 




— U.o 


4-Q9 1 


33 


Tempi. 


J 


754 


81±0 


05 


14 


395±0 


050 




— u.y 




27 


Tempi. 


11 


754 


47±0 


06 


14 


533±0 


050 




—0.7 


+92.0 


28 


Tempi. 


K 


755 


63±0 


05 


14 


532±0 


058 


SN 2006X 


-1.8 


+90.9 


8 


Tempi. 


B 


786 


15±0 


12 


15 


321±0 


003 


0.579±0.007 


-10.4 


+119.7 


41 


Spline 


V 


788 


91±0 


07 


14 


090±0 


006 


0.938±0.009 


-13.2 


+116.9 


41 


Spline 


u 


784 


55±0 


21 


17 


373±0 


013 


1.113±0.028 


-8.8 


+56.9 


23 


Spline 


9 


787 


79±0 


07 


14 


771±0 


003 


0.805±0.007 


-12.1 


+110.2 


39 


Spline 


r 


788 


77±0 


16 


13 


643±0 


006 


1.263±0.008 


-13.0 


+109.2 


39 


Spline 


i 


784 


54±0 


18 


13 


836±0 


005 


0.799±0.011 


-8.8 


+113.4 


39 


Spline 


Y 


777 


62±2 


06 


13 


346±0 


059 




-0.8 


+90.5 


42 


Spline 


J 


782 


77±0 


15 


12 


824±0 


009 


1.557±0.015 


-5.9 


+85.4 


43 


Spline 


11 


781 


69±0 


27 


12 


995±0 


007 


0.444±0.013 


-4.8 


+86.5 


42 


Spline 


K 


783 


95±0 


09 


12 


790±0 


031 


SN 2006ax 


+0.9 


+63.1 


18 


Tempi. 


B 


827 


40±0 


13 


15 


202±0 


003 


0.493±0.010 


-11.6 


+70.0 


29 


Spline 


V 


828 


51±0 


09 


15 


213±0 


002 


0.927±0.006 


-12.7 


+68.9 


26 


Spline 


u 


825 


70±0 


10 


15 


545±0 


007 


1.346±0.018 


-9.9 


+45.1 


21 


Spline 


9 


827 


65±0 


06 


15 


147±0 


003 


0.789±0.007 


-11.8 


+69.7 


26 


Spline 


r 


828 


74±0 


12 


15 


314±0 


003 


1.137±0.007 


-12.9 


+68.6 


26 


Spline 


i 


824 


06±0 


28 


15 


822±0 


006 


0.745±0.010 


-8.3 


+73.2 


26 


Spline 
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Table 2 — Continued 





^max 


m max 


Am t 


Epoch 




Fit 


Band 


-2,453,000 


(mag) 


(mag) 


Range 


Npts. 


Method 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 



Y 


823.46±0.17 


15 


841+0 


008 


0.043±0.017 


-5 


6 


+45 


4 


36 


Spline 


J 


823.85±0.19 


15 


756±0 


Oil 


1.494±0.030 


-6 





+45 





33 


Spline 


H 


823.58±0.13 


15 


914±0 


019 


0.410±0.045 


-5 


7 


+45 


3 


25 


Spline 


K 


824.53±0.07 


15 


862+0 


064 




+9 


3 


+19 


5 


6 


Tempi. 










SN 2006bh 














B 


833.73±0.06 


14 


456+0 


005 


0.729±0.007 


-4 


8 


+59 


4 


25 


Spline 


V 


834.88±0.08 


14 


461±0 


002 


1.154±0.007 


-5 


9 


+58 


3 


26 


Spline 


u 


831.88±0.08 


14 


867+0 


005 


1.614±0.010 


-2 





+61 


3 


21 


Spline 


9 


834.48±0.25 


14 


413±0 


006 


1.143±0.007 


-5 


5 


+58 


7 


24 


Spline 


r 


835.03±0.11 


14 


497+0 


006 


1.676±0.008 


-6 


1 


+58 


1 


26 


Spline 


i 


830.36±0.11 


14 


975±0 


004 


1.322±0.009 


-1 


4 


+62 


8 


25 


Spline 


Y 


830.12±0.04 


14 


881±0 


050 




-0 


2 


+39 


2 


24 


Tempi. 


J 


830.56±0.04 


14 


801±0 


050 




-0 


7 


+38 


7 


21 


Tempi. 


H 


830.26±0.04 


14 


977±0 


050 


SN 2006eq 


-0 


3 


+39 


1 


19 


Tempi. 


B 


977.11±0.35 


18 


593+0 


023 




+4 


8 


+49 


9 


14 


Tempi. 


V 


979.02±0.35 


18 


328+0 


013 




+2 


8 


+53 


1 


16 


Tempi. 


u 


975.45±0.35 


19 


339+0 


142 




+6 


5 


+29 


5 


6 


Tempi. 


9 


977.80±0.35 


18 


451±0 


014 




+4 





+54 


4 


17 


Tempi. 


r 


979.08±0.35 


18 


241±0 


009 




+2 


7 


+66 


8 


18 


Tempi. 


i 


976.88±0.35 


18 


653±0 


014 




+5 





+69 


1 


18 


Tempi. 


Y 


975.41±0.35 


18 


455±0 


025 




+5 


4 


+43 


2 


12 


Tempi. 


J 


976.21±0.35 


18 


459+0 


083 




+4 


6 


+42 


4 


9 


Tempi. 


H 


977.27±0.35 


18 


281±0 


138 


SN 2006gt 


+3 


5 


+41 


3 


9 


Tempi. 


B 


1004.01±0.28 


18 


343±0 


014 




-2 


2 


+59 


9 


12 


Spline 


V 


1005.53±0.37 


18 


112±0 


009 




-3 


7 


+58 


4 


13 


Spline 


u 


1001.53±0.12 


19 


048+0 


048 




+0 


1 


+18 





5 


Tempi. 


9 


1004.63+0.18 


18 


183±0 


009 




-2 


8 


+59 


3 


12 


Spline 


r 


1005.65±0.33 


18 


153±0 


Oil 




-3 


8 


+58 


3 


13 


Spline 


i 


1002.49±2.22 


18 


464+0 


031 




-0 


8 


+61 


4 


13 


Spline 


Y 


1002.92+0.12 


18 


387±0 


024 




-0 


3 


+47 


9 


20 


Tempi. 


J 


1003.29±0.12 


18 


358+0 


053 




-0 


6 


+39 


1 


16 


Tempi. 


H 


1004.08±0.12 


18 


236+0 


088 




-1 


4 


+38 


3 


11 


Tempi. 










SN 2006mr 














B 


1050.94±0.12 


15 


402+0 


014 


1.425±0.029 


-4 


2 


+96 


1 


30 


Spline 


V 


1052.83+0.08 


14 


612±0 


010 


1.761±0.019 


-6 


1 


+94 


2 


31 


Spline 


u 


1049.46±0.18 


16 


616±0 


017 


2.068±0.054 


-2 


8 


+50 


9 


19 


Spline 


9 


1051.37±0.11 


15 


020±0 


011 


1.647±0.017 


-4 


7 


+95 


6 


31 


Spline 


r 


1053.72±0.07 


14 


378±0 


014 


2.407±0.019 


-7 





+93 


3 


32 


Spline 


i 


1054.59±0.17 


14 


503±0 


009 


2.328±0.028 


-7 


9 


+92 


4 


31 


Spline 


Y 


1056.45+0.19 


14 


014±0 


007 


2.404±0.022 


-8 


7 


+80 


7 


29 


Spline 


J 


1054.12±0.12 


14 


056±0 


007 


3.486±0.119 


-6 


4 


+83 





19 


Spline 


H 


1055.31±0.07 


13 


850±0 


006 


3.299±0.081 


-7 


6 


+81 


8 


27 


Spline 
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Table 2 — Continued 





^max 


tn™ 


Amt 


Epoch 




Fit 


Band 


-2,453,000 


(mag) 


(mag) 


Range 


Npts. 


Method 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 






SN 2006py 








B 


1070.57+0.32 


18.140±0.016 




-1.1 +22.2 


7 


Tempi. 


V 


1072.16±0.32 


18.019±0.013 




-2.8 +20.6 


7 


Tempi. 


u 


1069.35±0.32 


18.594±0.052 




+0.2 +19.3 


5 


Tempi. 


9 


1071.16±0.32 


18.050±0.010 




-1.7 +17.4 


6 


Tempi. 


r 


1072.56±0.32 


18.058±0.009 




-3.2 +20.1 


7 


Tempi. 


i 


1068.11±0.32 


18.556±0.018 




+1.5 +24.9 


7 


Tempi. 



Note. — SNe with names in bold face belong to the "best-observed" group (see text). 
Columns: (1) Filter name; (2) Date of maximum; (3) Apparent peak magnitude (K- 
corrected); (4) Decline rate measured as the difference in magnitudes at an epoch t and 
maximum light in each band (t = 10 d for B; t = 15 d for ug; t = 20 d for V; t = 30 d 
for r; t = 35 d for iY; t = 40 for JH); (5) Epochs covered by data in rest-frame days 
with respect to the time of maximum; (6) Number of data points; (7) Fitting method 
(spline or template; see text) 
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Table 3. Pseudocolors versus Ami5(B) relation parameters 



Pseudocolor 


a 


b 


rms 


iVsNc 


(1) 


(2) 


(3) 


(4) 


(5) 


^max ymax 


0.30±0.04 


0.59±0.25 


0.14 


13 


^max ymax 


-0.02±0.01 


0.12±0.05 


0.06 


14 


^max _ ^,max 


-0.14±0.02 


0.24±0.06 


0.06 


14 


^max ^.max 


-0.09±0.02 


0.29±0.07 


0.07 


14 


ymax ^max 


-0.61±0.02 


0.27±0.08 


0.09 


14 


^max _ ymax 


-0.67±0.03 


0.51±0.17 


0.14 


13 


ymax jmax 


-0.63±0.02 


0.55±0.10 


0.09 


14 


ymax j^"max 


-0.77±0.03 


0.54±0.11 


0.14 


14 


ymax j^~max 


-0.59±0.05 


0.14±0.18 


0.14 


8 



Note. — Fits of the type: Pseudocolor = a + b [Amis(i?) — 
1.1] valid over the range 0.85 < Ami5(i?) < 1.76 mag. 
Columns: (1) Pseudocolor name; (2) Fit intercept; (3) Fit 
slope; (4) rms of fit in magnitudes; (5) Number of SNe used 
in fit. 



Table 4. Host-galaxy color excesses from maximum light colors 



SK 


E(B - V) 


E(B - V) 


E(u - V) 


E{g - r) 


E(B - r) 


E{V - i) 


E{V - Y) 


E(V - J) 


E(V - H) 


E{V - K) 




Tail 


Max 


Max 


Max 


Max 


Max 


Max 


Max 


Max 


Max 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


0) 


(10) 


(11) 


2004ef 


0.125(017) 


0.089(006) 


0.343(014) 


0.117(005) 


0.137(007) 


-0.037(006) 


0.117(034) 


0.024(080) 




-0.176(473) 


2004oo 


0.082(018) 


0.020(016) 


0.070(026) 


-0.023(018) 


0.006(015) 


-0.004(018) 


-0.244(021) 


0.054(021) 


-0.093(022) 


0.013(017) 


2004ey 


0.039(005) 


-0.061(008) 


0.006(015) 


-0.004(006) 


-0.043(008) 


-0.054(008) 


-0.111(020) 


0.089(051) 


-0.055(024) 


-0.159(046) 


2004gc 


0.200(021) 


0.149(037) 




0.206(038) 


0.235(041) 


0.078(038) 










2004gs 


0.080(010) 


0.121(007) 


0.300(020) 


0.092(007) 


0.138(008) 


0.080(009) 




0.063(019) 






2004gu 


0.257(017) 


0.200(024) 


0.549(045) 


0.263(013) 


0.358(025) 


0.128(019) 










2005A 




1.005(011) 


2.309(047) 


1.078(013) 


1.451(014) 


0.611(018) 


1.414(031) 


1.540(031) 


1.483(039) 




2005M 


0.091(007) 


0.050(005) 


0.144(019) 


0.028(008) 


0.067(008) 


0.082(006) 


0.137(015) 


0.210(048) 


0.115(020) 


0.144(036) 


2005W 




0.352(011) 


0.326(021) 


0.182(015) 


0.428(013) 


0.114(015) 










2005ag 


0.064(026) 


0.050(009) 




0.076(007) 


0.069(008) 


0.227(021) 


-0.065(022) 


-0.041(038) 


0.272(056) 




2005al 


-0.050(004) 


-0.045(009) 


-0.008(025) 


-0.037(012) 


-0.056(014) 


-0.135(014) 


-0.101(023) 


-0.020(025) 


-0.174(026) 


-0.227(064) 


2005am 


0.050(005) 


0.014(005) 


0.034(011) 


0.027(007) 


0.023(007) 


-0.046(011) 




0.044(017) 


0.046(014) 


0.007(043) 


2005be 




-0.012(047) 


-0.186(069) 


-0.091(047) 


-0.092(050) 


-0.183(042) 










2005bg 


-0.005(028) 


0.063(062) 


0.152(074) 


0.040(066) 


-0.000(066) 


-0.137(071) 










2005bo 


0.132(110) 


0.277(020) 


0.455(060) 


0.286(023) 


0.360(026) 


0.232(038) 










2005el 


-0.050(008) 


-0.064(010) 


-0.168(011) 


-0.153(014) 


-0.121(013) 


-0.044(009) 


0.071(013) 


-0.085(012) 


-0.039(013) 


-0.096(135) 


2005eq 


0.099(012) 


0.126(009) 


0.263(018) 


0.163(009) 


0.205(010) 


0.054(051) 


0.267(051) 


0.247(031) 


0.142(041) 


0.274(073) 


2005hc 


0.078(017) 


0.035(007) 


0.245(013) 


0.059(008) 


0.070(007) 


-0.105(011) 


0.333(051) 


0.140(047) 


0.205(103) 




2005iq 


0.028(023) 


-0.033(012) 


-0.063(016) 


-0.069(010) 


-0.072(011) 


0.001(015) 


0.204(051) 


-0.030(047) 


0.142(051) 




2005ir 


-0.010(165) 


0.097(038) 


0.209(087) 


0.135(036) 


0.107(044) 


-0.110(066) 










2005kc 




0.172(008) 


0.422(014) 


0.239(008) 


0.295(008) 


0.273(009) 


0.482(013) 


0.552(013) 


0.515(040) 




2005ke 


0.036(005) 




















2005ki 


-0.048(005) 


-0.003(007) 


-0.074(015) 


-0.043(021) 


-0.031(021) 


-0.006(017) 


-0.037(017) 


-0.030(073) 


0.027(030) 




20051u 


0.217(021) 


0.234(023) 


0.544(063) 


0.257(020) 


0.367(023) 


0.147(022) 










2005na 


0.113(011) 


-0.017(011) 


-0.152(074) 


0.022(010) 


0.004(013) 


0.071(051) 










2006D 


0.014(004) 


0.074(007) 


0.252(015) 


0.054(008) 


0.074(009) 


0.060(007) 


-0.051(050) 


0.134(051) 


0.122(051) 


0.063(058) 


2006X 




1.221(007) 


2.943(017) 


1.242(008) 


1.737(008) 


0.838(009) 


1.365(060) 


1.836(013) 


1.796(012) 


1.815(032) 


2006ax 


-0.011(005) 


-0.032(004) 


0.013(011) 


-0.055(006) 


-0.058(006) 


-0.023(008) 


-0.003(011) 


0.028(014) 


-0.007(021) 


-0.188(064) 


2006bh 


-0.009(004) 


-0.052(006) 


-0.122(009) 


-0.048(009) 


-0.074(009) 


-0.017(006) 


0.039(050) 


0.055(050) 


0.012(050) 




2006eq 




0.168(027) 


0.301(146) 


0.162(021) 


0.217(029) 


0.080(024) 


0.163(040) 


0.076(089) 


0.383(142) 




2006gt 


-0.021(063) 


0.144(017) 


0.243(053) 


-0.006(016) 


0.070(020) 


0.065(034) 


0.036(031) 


-0.014(057) 


0.239(090) 




2006mr 


-0.025(013) 




















2006py 




0.085(021) 


0.219(056) 


0.083(015) 


0.109(020) 


0.026(024) 











Note. — Columns: (1) SN name; (2) Host-galaxy reddening measured from the (B — V) tail (§ \3.2\ ; (3)-(ll) Host-galaxy reddening measured from pseudocolor at maximum 

light (§E3J. 
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Table 5. Ry values from fits to E(V — X) max versus E(B — V) 





E(V-i) 


E(V - Y) 


E{V -J) 


E(V - H) 


Weighted 


Sample 


Max 


Max 


Max 


Max 


Average 


Whole 


1.1 ±0.2 


1.6 ±0.2 


1.8 ±0.1 


1.7±0.1 


1.7 ± 0.1 


Excl. 2005A and 2006X 


3.8 ±1.5 


3.5 ± 1.4 


3.1 ±0.7 


3.1 ±0.7 


3.2 ±0.4 
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Table 6. Conversion of observed to true color excesses 



Color Excess 


Correction Factor Cx 


E(V - u) 


1.0576 


E(V - B) 


1.1023 


E(V - g) 


1.1855 


E(V - r) 


1.2028 


E(V-i) 


1.0218 


E(V - Y) 


1.0370 


E(V - J) 


1.0142 


E(V - H) 


0.9951 


E(V - K„) 


0.9818 



Note. — E(V - X) tIUC = C x * E(V - 

observed 



Table 7. Adopted distances for three SNe with z < 0.01 





Host 




Distance 




SN 


Galaxy 


(mag) 


Method 


Reference 


2005ke 


NGC 1371 


31.84 ±0.08 


SBF 


Tonrv et al. (2001) a 


2006X 


NGC 4321 


30.91 ±0.14 


Cephcids 


Freedman et al. (20011 


2006mr 


NGC 1316 


31.59 ±0.08 


SBF 


Cantiello et al. (2007) 



a NGC 1371 is a member of the Eridanus group/cluster. iTonrv et al.l (|200lh 
give a distance modulus of 32.00 ± 0.08 mag for this group based on SBF 
measurements of 7 members. Subtrac ting 0.16 mag ( Jensen et al. 2003f) to 
put this on the lFreedman et al.l (|2001r ) scale gives 31.84 ± 0.08 mag. 



Table 8. Fits of peak magnitudes versus Am\^{B) and pseudocolor 



Fit 


Filter 


Pseudocolor 








Rv 


""SAT 


RMS 






No. 


X 


(Y-Z) 


M x (0) 


bx 




(CCM+O) 


[mag] 


[mag] 


N S Ne 


Sample 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


1 


B 


(B-V) 


-19.10±0.02 


0.64±0.10 


2.81±0.11 


1.52±0.10 


0.14 


0.17 


31 


whole 


2 


B 


(B-V) 


-19.09±0.02 


0.74±0.10 


2.74±0.11 


1.46±0.09 


0.12 


0.15 


26 


best-observed 


3 


B 


(B-V) 


-19.10±0.01 


0.61±0.11 


2.79±0.12 


1.50±0.10 


0.12 


0.16 


28 


excluding fast-decliners 


4 


B 


(B-V) 


-19.09±0.01 


0.74±0.11 


2.74±0.11 


1.46±0.10 


0.09 


0.14 


23 


best-observed excl. fast-decliners 


5 


B 


(B-V) 


-19.07±0.02 


0.81±0.12 


2.24±0.37 


1.03±0.32 


0.09 


0.14 


21 


best-observed and (B — V) < 0.4 


6 


J 


(V-J) 


~18.33±0.03 


0.65±0.10 


0.20±0.07 


1.60±0.32 


0.08 


0.14 


21 


best-observed 


7 


J 


(V-J) 


-18.36±0.03 


0.52±0.09 


0.13±0.06 


1.32±0.22 


0.04 


0.12 


19 


best-observed excl. fast-decliners 


8 


J 


(V-J) 


-18.44±0.14 


0.57±0.14 


0.01±0.24 


0.96±0.59 


0.05 


0.12 


17 


best-observed and (V — J) < 0.0 



Note. — Fits of the type: nx = m™ ax - M x (0) - b x [Am 15 (B) - 1.1] - /3^ z (Y-Z). 
Columns: (1) Fit identifier; (2) Filter corresponding to m^ ax ; (3) Color; (4) Absolute magnitude for Ami5(B) = 1.1 and zero (Y — Z) color; (5) Luminosity- 
decline rate slope; (6) Luminosity-color slope; (7) Corresponding parameter Ry of the CCM+O reddening law; (8) Resulting intrinsic dispersion of SN 
data; (8) RMS of fit in magnitudes; (10) Number of SNe used in fit; (11) Sample of SNe used in fit (see text). 



Table 9. Fits of peak magnitudes versus Amis (B) and reddening 



r It 


Filter 


Color excess 








Hy 


o"5Af 


Jxlvlo 






INO. 


A 


E(Y - 


Z) 


M x {v) 


I, 

ox 


r/YZ 

K x 




[mag] 


[mag] 


AT 


Sample 


(i) 


( 2 ) 


(3) 




( 4 ) 


(5) 


(6) 


(<) 


(8) 


(9) 


(10) 


I 11 ) 


i 


B 


E(B - 


V) 


-19.07±0.01 


1.03±0.25 






0.19 


0.19 


12 


low reddening 


2 


B 


E(B - 


V) 


-19.39±0.02 


0.97±0.41 


4.53 


3.10 


0.54 


0.53 


28 


whole 


3 


B 


TP 1 D 
tj\D — 


V) 


-19.15±0.02 


0.94±0.11 


2.79±0.12 


1.50±0.11 


0.13 


0.16 


28 


whole 


4 


B 


E(B - 


V) 


-19.14±0.01 


1.07±0.11 


2.74±0.11 


1.46±0.10 


0.10 


0.14 


23 


best-observed 


5 


B 


E{B - 


V) 


-19.14±0.03 


0.94±0.12 


2.75±0.38 


1.46±0.33 


0.13 


0.16 


26 


excluding 05A and 06X 


6 


B 


E(B - 


V) 


-19.11±0.02 


1.08±0.11 


2.22±0.37 


1.01±0.31 


0.09 


0.14 


21 


best-observed excl. 05A and 06X 


7 


V 


E(B - 


V) 


-19.12±0.01 


0.95±0.11 


1.74±0.11 


1.46±0.10 


0.09 


0.14 


23 


best-observed 


8 


u 


E(u - 


V) 


-18.83±0.02 


1.59±0.12 


1.74±0.05 


1.10±0.07 


0.10 


0.13 


22 


best-observed 


9 


9 


E(g- 


r) 


-19.17±0.01 


1.32±0.12 


2.21±0.11 


1.49±0.12 


0.11 


0.14 


23 


best-observed 


10 


r 


E(B - 


■r) 


-19.03±0.01 


1.02±0.11 


0.89±0.08 


1.44±0.11 


0.10 


0.13 


23 


best-observed 


11 


i 


E(B - 


V) 


-18.50±0.02 


0.71±0.13 


1.08±0.12 


1.71±0.13 


0.12 


0.15 


23 


best-observed 


12 


Y 


E(B - 


V) 


-18.46±0.01 


0.44±0.14 


0.54±0.12 


1.78±0.19 


0.10 


0.15 


18 


best-observed 


13 


J 


E(V - 


J) 


-18.44±0.01 


0.58±0.09 


0.14±0.06 


1.37±0.23 


0.03 


0.12 


18 


best-observed 


14 


H 


E(y - 


H) 


-18.38±0.02 


0.32±0.18 


0.20±0.09 


2.72±1.79 


0.12 


0.16 


16 


best-observed 


15 


K 


E(V - 


K) 


-18.43±0.03 


0.81±0.29 


0.15±0.10 


4.11±8.32 


0.10 


0.17 


9 


best-observed 



Note. — Fits of the type: = ~ M x(0) ~ b x [Am 15 (B) - 1.1] - R Y X Z E(Y - Z). 

Columns: (1) Fit identifier; (2) Filter corresponding to mj£ ax ; (3) Color excess; (4) Absolute magnitude for Amis(B) = 1.1 and no reddening; (5) 
Luminosity-decline rate slope; (6) Total-to-selective absorption coefficient R x z = Ax/E(Y — Z); (7) Corresponding parameter Ry of the CCM+O 
reddening law; (8) Resulting intrinsic dispersion of SN data; (9) RMS of fit in magnitudes; (10) Number of SNe used in fit; (11) Sample of SNe used in fit 
(see text). 



Table 10. Precision of the SNe la as standard candles in the NIR 



Band 
X 


M max 

(avg.) 


Adopted 

r>YZ 

H x 


Rv 
(CCM+O) 


RMS 
[mag] 


N S Ne 


Sample 




Y 
Y 
Y 


-18.44 ±0.07 
-18.47 ±0.08 
-18.43 ±0.08 


R\ Y = 0.32 
R\- Y = 0.60 
R\ Y = 0.0 


1.7 
3.1 
0.0 


0.24 
0.26 
0.19 


17 
15 
15 


best-observed 
excluding 05A and 
excluding 05A and 


06X 
06X 


J 
J 
J 


-18.43 ±0.07 
-18.44 ±0.07 
-18.42 ±0.07 


R V j J = 0.13 
Ry = 0.39 

Ry = o.o 


1.3 
3.1 
0.0 


0.18 
0.20 
0.18 


18 
16 
16 


best-observed 
excluding 05A and 
excluding 05A and 


06X 
06X 


H 
H 
H 


-18.42 ±0.08 
-18.43 ±0.08 
-18.40 ±0.08 


R V H H = 0.19 
R V H H = 0.22 
R V H H = 0.0 


2.4 
3.1 
0.0 


0.19 
0.21 
0.19 


16 
14 
14 


best-observed 
excluding 05A and 
excluding 05A and 


06X 
06X 


K s 
K s 
K s 


-18.47 ±0.12 
-18.49 ±0.13 
-18.49 ±0.13 


R V K K = 0.12 
rvk = 0.13 
R V K K = 0.0 


2.7 
3.1 
0.0 


0.27 
0.29 
0.27 


9 
8 
8 


best-observed 
excluding 06X 
excluding 06X 
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Table 11. Comparison of Absolute Magnitudes in JHK S 



Band 


Maximum 
Krisciunas et al. (2004b) 


This Paper 


^{^max ) 

Wood-Vasev et al. (2008) 


This Paper 


J 


-18.61 ±0.03 


-18.43 ±0.07 


-18.29 ±0.09 


-18.42 ±0.02 


H 


-18.28 ±0.03 


-18.42 ±0.08 


-18.08 ±0.09 


-18.23 ±0.03 


K s 


-18.44 ±0.03 


-18.47 ±0.12 


-18.30 ±0.11 


-18.30 ±0.05 
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Table 12. Measurements of Ry for Individual SNe la 



SN 


Am 15 (B) 


E(B - V) 


Rv 


Reference 


1999cl 


1.29 ±0.08 


1.24 ±0.07 


1.55 ±0.08 


Krisciunas et al. (20061 


2001el 


1.15 ±0.04 


0.21 ±0.05 


2.15 ±0.23 


Krisciunas et al. (20071 


2002cv 


1.46 ±0.17 


5.45 ±0.28 


1.59 ±0.07 


Elias-Rosa et al. (2008) 


2003cg 


1.25 ±0.05 


1.33 ±0.11 


1.80 ±0.19 


Elias-Rosa et al. (2006) 


2005A 


1.34 ±0.03 


1.11 ±0.07 


1.68 ±0.10 


This paper 


2006X 


1.31 ±0.05 


1.42 ±0.04 


1.48 ±0.06 


Wang et al. (2008a) 




1.29 ±0.05 


1.34 ±0.07 


1.55 ±0.07 


This paper 



Note. — The value s Ami^(B) given her e have been corrected for reddening 
using equation (6) of Phillips et al. ( 1999h . 
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Table 13. Reddening-law coefficients for a SN la at maximum 



CCM CCM+O CCM+O (Obs.) 



Filter 


ax 


bx 


ax 


bx 


ax 


bx 


u 


0.990 


1.605 


1.006 


1.480 


0.959 


1.849 


B 


1.016 


0.860 


1.026 


0.799 


1 


1 


9 


1.017 


0.540 


1.020 


0.525 


1.003 


0.657 


V 


1 





1 





1 





r 


0.939 


-0.245 


0.942 


-0.284 


0.952 


-0.356 


i 


0.807 


-0.499 


0.797 


-0.449 


0.813 


-0.564 


Y 


0.539 


-0.502 


0.538 


-0.503 


0.557 


-0.636 


J 


0.398 


-0.370 


0.398 


-0.371 


0.412 


-0.469 


H 


0.259 


-0.243 


0.259 


-0.244 


0.268 


-0.308 


K 


0.166 


-0.153 


0.165 


-0.153 


0.171 


-0.193 
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50 100 50 100 50 100 

rest-frame days since B maximum 

Fig. 1. — Spline fits to the ugriBVY J H light curves of the SNe used as templates. The light 
curves are sorted by Ami^(B) and are plotted relative to the magnitude at maximum light in each 
band, shifted by 0.5 mag with respect to each other. The values of Ami^(B) label each curve at 
the right end. 
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Fig. 2. — Comparison of decline rates Ami5(S) with similar quantities obtained from other optical 
bands and epochs since maximum light. Decline rates Amt(X) are measured from the light curve 
in X as the difference in magnitudes between t rest-frame days since maximum light (in that band) 
and the peak magnitude. Only results of spline-function fits to SNe with pre-maximum coverage 
are shown for each band. Except where plotted, error bars are smaller than the points. 



-59- 



1.5 - 

1 r 

2 1 
Si. 5 

lO 
■<—< 

B 

< 1 

2 
1.5 
1 



1 I , I 


1 1 1 1 1 1 1 1 1 1 


i i i i i 


: Am 35 (Y) 


• • 




- 




- 

_ 


- 

- • 




- 


• 




- 


— — 1 — 1 — h ~ 
: Am 40 (J) 


H — — 1 — 1 — 1 — 1 — — 1 — 1 — 1— 


H — — 1 — 1 — h 


- 


-•- • 


- 


- 


• 


- 




• 




- 


• 

H — — 1 — 1 — 1 — 1 — — 1 — 1 — h- 


- 

H — — 1 — 1 — h 


— — 1 — 1 — h_ 
: Am 40 (H) 


O 


• 




«* 




• 






• 






i , , , 


i 1 i i i i 1 i i i 


i 1 i i i 





1 2 


3 



Am t (X) 

Fig. 3. — Comparison of decline rates Amis (5) with similar quantities obtained from YJH bands. 
Decline rates Amt(X) are measured from the light curve in X as the difference in magnitudes 
between t rest-frame days since maximum light (in that band) and the peak magnitude. Only 
results of spline- function fits to SNe with pre-maximum coverage are shown for each band. Except 
where plotted, error bars are smaller than the points. 
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Fig. 4. — Differences in rest-frame days between the time of maximum in ugriVY JHK S with 
respect to B. Only results of spline-function fits to SNe with pre-maximum coverage are shown for 
each band. 
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Fig. 5. — (left) Plot of the difference in the peak magnitudes derived from template versus spline 
fits as a function of filter for the subset of best-observed SNe. (right) Difference in Ami^(B) derived 
from template versus spline fits for the same sample. 
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Fig. 6. — Timing and strength of the secondary maximum in i as a function of the decline-rate 
parameter /\m\^{B) for the sample of best-observed SNe with pre-maximum data in each band. 
The top panel shows the phase of the secondary maximum with respect to the £>-band maximum. 
The middle panel displays the difference in magnitudes between the secondary maximum and the 
local minimum between the primary and secondary maxima. The bottom panel shows the difference 
in magnitude between the primary and secondary maxima. 
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0.2 0.4 0.6 

(i)20-40 

Fig. 7. — The mean i-band flux 20 to 40 days after the time of i?-band maximum plotted versus 
Am\5(B) for the best-observed subsample of CSP SNe. Solid symbols correspond to SNe with peak 
magnitudes measured from spline fits; open symbols are SNe with peak magnitudes derived from 
template fits. The solid line corresponds to the relation for the SNOOPy templates. Four SNe with 
similar Amis(-B) but discrepant average i-band fluxes at the second maximum are labeled. 
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Fig. 8.— (left) Comparison of the B, V, and i light curves of SNe 2006D, 2006bh, 2004eo, and 
2004ef. The data for SNe 2006D, 2004eo, and 2004ef have been shifted to coincide at maximum 
with the observations of SN 2006bh. All four SNe have very similar decline rates in the range 
Ami^(B) =1.37-1.42 mag. Note the close similarity of the B and V light curves, but the real 
differences in the strength and morphology of the secondary maxima in the i-band. (right) The 
difference between the maximum-light magnitudes in the BVi bands as measured from template 
and spline fits for SNe 2006D, 2006bh, 2004eo, and 2004ef, plotted as a function of the mean z-band 
flux 20 to 40 days after the time of .B-band maximum. 
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Fig. 9. — Late-time B — V color evolution of SNe assumed to have low or no host-galaxy reddening. 
Observed colors were corrected for Galactic reddening. The solid line shows a linear fit to the data 
in the range 30 < ty < 80 days (see text). The rms scatter about this fit is 0.077 mag. 
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Fig. 10. — Equivalent- width measurements (or upper limits) of the Na I D lines plotted versus the 
color excess in B — V as measured from the lLiral (|1995l ) law. The rms dispersion of 0.077 mag in the 
Lira-law fit (equation [2]) has been added in quadrature to the errors in E(B — V)taii- Upper limits 
to the equivalent width are indicated as "U." The red symbols with dashed error bars correspond 
to the low-reddening sample, while the blue symbols with solid error bars show all the other SNe in 
our sample with measurements of E(B — V)tai\ an d the Na I D equivalent width (or upper limit). 
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Fig. 11. — Pseudocolors at maximum light corrected for Galactic reddening, as a function of 
Ami5(£?), for several filter combinations. Filled circles correspond to the low-reddening subsample 
of SNe (see § I3.1j) . The dashed lines are linear fits to the data of the low-reddening SNe in the 
range of 0.8 < Ami5(B) < 1.7 mag. 
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Fig. 12. — Comparison of color-excess estimates for B — V from measurements at maximum light 
[E(B — y) max ] and at late epochs [E(B — V) ta ii]- The difference between these two quantities is 
plotted vs. E(B — F) max , B max — V ma _ x , Amxs(5), and the absolute magnitude in B corrected for 
decline rate (using Fit 1 of Table EJ see § I4.2[) . The dotted lines in the top two panels indicate the 
weak correlations that would be expected due to the fac t that the color exce sses given in this paper 
are "observed" [as opposed to "true" color excesses (see iPhillips et al.lll999l )]. and are measured at 
very different epochs in the color evolution of the SN. 
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Fig. 13. — Comparison of color-excess estimates E(V — X) max , for X = iYJH, with E(B — V) ma , x 
for the best-observed SNe. The solid red lines represent the slope which corresponds to the average 
fit value of Ry = 1-7, found using the whole set of points. The dashed blue lines indicate the slope 
predicted by a value of Ry = 3.2, which is the averaged fit value found when excluding the two 
highly reddened SNe (the two points farthest to the right in the plots.) 



-70- 




Fig. 14.— Color excesses E(V - X\), for bands X x = ugriBY JHK S , for SN 2005A (left panel) 
and SN 2006X (right panel). The solid lines show the best-fit CCM+O laws, whereas the dotted 
lines sho w the CCM+O model for Ry = 3.1. The dashed lines correspond to fits of the power-law 
model by iGoobarl (|2008l ). valid for LMC-type dust. The latter provides a substantially better fit of 
the u-band reddening for both SNe, as compared to the CCM+O model. 
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Fig. 15.— Color excesses E(V - X x ), for bands X x = ugriBYJHK s , for SN 2005kc (left panel) 
and SN 2006eq (right panel). See the Figure RBI caption for further details. 
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Fig. 16. — Sample fits of absolute peak magnitude versus decline rate (Amis(i?)) and color for SNe 
with reliable distances. One fit per row is displayed. The left panels show the relations between 
color-corrected absolute peak magnitude and decline rate. The right panels show the absolute peak 
magnitude vs. color relations, corrected for decline rate. Red circles denote the fast-declining SNe; 
red triangles mark the very red SNe 2005A and 2006X. The resulting fits are shown with solid lines. 
The dotted lines in these panels show the relations predicted by reddening with the CCM+O law 
and Ry = 3.1. The number of SNe used in the fit, the rms scatter, the intrinsic dispersions (<tsn)j 
and the predicted Ry are labeled. Rows : (a) B vs. Amis (B) and B — V for the whole sample; (b) 
B vs. Am\5(B) and B — V for the subsample of best-observed SNe; (c) J vs. Amis (B) and V — J 
for the subsample of best-observed SNe. 
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Fig. 17. — Fits of absolute peak magnitude versus decline rate (Ami5(B)) and reddening (measured 
as a color excess). One fit per row is displayed. The left panels show the relations between 
reddening-corrected absolute peak magnitude and decline rate. The right panels show the absolute 
peak magnitude vs. color excess relations, corrected for decline rate. The symbols for points and 
lines, and the labels, are equivalent to those used in Figure [TBI For comparison, we ad ded to the left 



panels of the upper two rows a dashed line indicating the quadratic relation fits by iPhillips et al 



( 19991 ) . shifted to the distance scale adopted in this work which corresponds to Hq = 72 km s _1 
Mpc . Note that fast-declining SNe are not included in these plots. Rows: (a) B vs. Amis(5) 
and E(B — V) for the whole sample; (b) B vs. Amis(-B) and E(B — V) for the subsample of 
best-observed SNe; (c) J vs. Amis(-B) and E(V — J) for the subsample of best-observed SNe. 
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Fig. 18. — Absolute iYJH light curves for the subsample of best-observed SNe. The observed 
magnitudes are corrected for K-corrections, extinction in the Galaxy (assuming Ry = 3.1), and 
extinction in the host galaxies (using the FCC values of Fits. 11 to 14). The times of observation 
are put in the rest frame of the SNe, using the heliocentric redshifts, and referred to the time of B 
maximum. The distances to the SNe are determined as described in § SJ The data from different 
SNe are shown with different symbols and the template light-curve fits are included for guidance. 
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Fig. 19. — Hubble diagram for the sample of best-observed SNe out to z ~ 0.08. The top panel 
shows the combined distance moduli for each SN from the fits of absolute peak magnitude vs. 
decline rate and reddening in ugriBVY J H . The solid line shows the adopted concordance model 
of equation (0). The dotted lines represent the spread predicted by a dispersion of 382 km s _1 in 
the measured redshifts due to peculiar motions of the host galaxies. The bottom panel shows the 
residuals with respect to the model fit. 
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Fig. 20. — Residuals in the distance moduli calculated in band X, where X = ugriYJHK, plotted 
versus the residuals in the 5-band distance moduli. Note the significant correlation between these. 



